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Mixed Path Ground Wave Propagation: 
1. Short Distances 


James R. Wait 


An expression is derived for the mutual impedance between two short vertical antennas 
on a flat earth with a straight boundary separating two media of differing electrical constants. 
After making some approximations that are valid at low and medium frequencies and where 
the antennas are not near the boundary, the integral formula for the field is evaluated for a 
wide range of the parameters. The numerical results computed in this paper are shown to be 
in reasonably good agreement with experiment. Finally, the effect of the obliqueness of the 
boundary is considered by a refinement of the stationary phase evaluation of the integrals. 


1. Introduction 


Considerable interest has been shown recently in the propagation of ground waves over an 
inhomogeneous conducting earth. For many purposes, such as estimating coverage of broad- 
cast transmitters, it is usually sufficient to assign an equivalent conductivity to the path. A 
simple and effective method for estimating this equivalent value from a conductivity profile of 
the path has recently been proposed by Suda [1].’ Another technique described by Kirke [2], 
which is similar, has been called the equivalent-distance method and is also simple to apply. 
The best known method, however, is due to Millington [3]. Although semiempirical it does 
appear to be valid for a wide range of frequencies and ground constants. It also predicts the 
recovery effect that occurs when the wave passes over a boundary from an area of poor 
conductivity to one of good conductivity. Millington does not provide a theoretical justification 
for his method, although he does indicate that his formulas for high frequencies are compatible 
with the expected behavior of the height-gain functions over the media on both sides of the 
boundary. Clemmow [4] in an elegant dissertation obtains a rigorous solution for a line source 
on flat earth parallel to a boundary separating two media. He makes a limited comparison of 
his formula with Millington’s method, and shows that the agreement is good. In a more 
recent theoretical approach to the subject, Bremmer [5] formulates the problem in terms of an 
integral equation, which he solves by operational methods. After considerable manipulation he 
succeeds in showing that his result is mathematically equivalent to that of Clemmow. 
Bremmer considers several limiting cases and also establishes the validity of Millington’s 
formulas at high frequencies. 

It is the purpose of the present paper to extend Bremmer’s result, with particular attention 
being paid to the phenomena at low and medium radiofrequencies where it does not seem 
possible to obtain convenient series expansions for the field by the Bremmer method. For 
the sake of completeness, the problem is reformulated in terms of mutual impedance between 
two antennas located on a flat earth, with a straight boundary separating two homogeneous 
media. In this case there is no doubt as to whether the reciprocity theorem is satisfied or not. 


The line joining the two antennas makes some angle with the boundary that should not be 
near zero. An integral equation is obtained for the field which is similar to one formulated 


by Feinberg [6] for propagation over a rough ground. In the present paper, the integral 
equation is solved by numerical means for a range of parameters that are appropriate for short 
distances at low and medium frequencies. It is intended to extend the calculations in a later 
paper to situations where the earth’s curvature must be considered; that is, for distances in 
miles greater than about 50\’*, where \ is the wavelength in meters. 


e the literature references at the end of this paper 
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2. Formulation 


The surface of the earth is considered to be flat and defined by 0 in a simple cartesian 
coordinate system (z,y,2). Points above the earth correspond to positive values of 2. The 
earth medium to the left of a boundary line defined by y=z2 tan @, as indicated in figure 1, has 
a conductivity o and dielectric constant «. The earth medium for points to the right of the 
boundary line has a conductivity of o; and dielectric constant ¢. Short hertzian dipole antennas 
are located at the points A and B, which can be located anywhere along the z axis just above 
the surface of the ground at z=0. 

It is now assumed that a current J, applied at the terminals of antenna A would produce 
electric and magnetic fields E., and H, for the case when the earth media were homogeneous 
with electrical constants o and «. <A time factor exp(i#t) is implied. The mutual impedance 
between antennas A and B for this homogeneous case is denoted as Z,).. When the earth becomes 
inhomogeneous in the manner described above, the fields of antenna A with the same current 
I, become F, and Hi, and the mutual impedance becomes Z;,. It is convenient to regard the 
change from Z,, to Z,) as a consequence of the changes of the currents within the earth when 
it is transformed from a homogeneous to an inhomogeneous state. These currents will result 
in electric and magnetic fields E.—E, and Hi—H, and a voltage J,(Z,,—Z») at the terminals 
of antenna B. A current J, is now considered to be applied at terminals of antenna B, and 
the resulting electric and magnetic fields are £, and H,, over the homogeneous earth of con- 
stants o and «. It then follows from Ballentine’s ‘corollary L’’ of the electromagnetic reci- 
procity theorem [7,8] that 


ats l ae — =< 
La—Lar LI, | |. (hx He— ha & Ay), dr dy, (1) 


where the integration extends over the whole ground plane S, and the subscript 2 indicates 
that the 2 or normal component of vector products is taken. 

A simplification is now made by introducing the concept of surface impedance [6,8], that 
is, the tangential electric and magnetic fields on the surface of the earth are assumed to be 
related by a complex constant of proportionality. More specifically, 


E.~—n1H, } 


where 7, the surface impedance, is assumed to be equal to the value obtained for a plane wave 
at grazing incidence on the flat earth. It is given by 


n 120 2/1 +8] ’ (3) 
a ar 


where B=22/wavelength, and a=(icuw—euw*)"*, with w=42%107-7. A subscript 1 is to be 
added to 7 when the electric constants are o,; and «. 


Equation (1) can now be written 


1,1,AZ=(n— m) | (H1i,-H,)drdy, (4) 
where AZ=Zi,—Z,» is the change of the mutual impedance from the situation of the homo- 
geneous earth of surface impedance y to the inhomogeneous model indicated in figure 1. The 
quantity /7,, is the tangential magnetic field of the antenna B over the homogeneous earth, 


and #/,, is the tangential magnetic field of the antenna A over the inhomogeneous earth. 








—, T> 





















0; 
(x,y, O) 











ao ane 
-_-_ 


“Boundary Line 


FIGURE 1. Schematic diagram of antennas A and B on a flat earth with a straight bound- 


ary separating two media with different conductivities. 


The integration now extends over the surface So to the right of the boundary line. Since 
H.,, is not actually known, eq (4) is actually a surface integral equation formulation of the 
problem. Before attempting to solve it, a simplification can be made if the principle of sta- 
tionary phase is utilized. 

Letting the hertzian antenna at B have an effective height A,, the tangential magnetic 
field at P, a distance r from B, considering the earth as homogeneous and flat, is given by 


ip Thy 


2ri 


H,, * (its 


=) Fir, 9), X1,), (5) 


where F(r,n) is Sommerfeld’s ground-wave attenuation factor defined by 


Fir, n) ] (2p P é w da, (6) 
e/ a=ip” 
where 
iBr B(, , B\ Brn 
pa Ts | 1+%3) _ 2 ( 302) 


is the numerical distance. Equation (6) is not exact, being an approximation valid for 


8<a’*. 7, and 2, are unit vectors in the directions of increasing 7 and z, respectively. The 
function F(r,n) has been treated numerically by Norton [9] and others. The tangential mag- 


netic field 14, at a distance PR from A is of the form 


IBIahe a l " > _ _ a 
deh © (A tagp ) FR, mm) inte), (7) 


where F’(2,y,m) is some function of PR, 7, and m, and can be expected to be slowly varying 


compared to e~*. 7, is a unit vector in the direction of increasing PR. 














Using eq (5) and (7) and denoting the angle between the vectors 7, and iz as 4, it follows that 
the expression for the mutual increment AZ can be written 
, Bhh, (Ce #°r* l | , 
AZ ( ) " || > 1+. ie : ) Fr, ) x I ‘(r, »m)) COS bdardy. R 
1—™) ae? JJ. rh ( > T aR u) 1, odrdy, (8) 
where 


r=((7o+2)’+y'*|" and 
R=|(Ry—xzr)?+ y’|"*. 
The major contribution to the integrand occurs when the phase of the exponential term is 


nearly constant, since the others are relatively slowly varying. With this in mind, the 
exponent is expanded in a powers series in ¥’ as follows: 


r+ Rao Rote ( = re ) +( i - )] for Ry >r >0 
r+ R=ro— Rot 21 +0 ( r r )+( r “E | spintiiee “) 


and 


where terms in y’, y’, etc., are neglected. A further approximation to the integrand in eq (8) 
is to replace (1+-1/7/8r) and (141/78) by unity. This will be justified if the antenna A or B 
is not near the boundary (i. e., By and Bre>1). This is essentially equivalent to stating 
that the effect of the induction field of the antennas is not considered so far as the boundary 
is concerned. 

The integral expression for AZ now has the following approximate form: 


. ( \B hgh, e P09") C8 F(rotaz,n)F’(Ro—x,0,m) abe 2 
AZ 4. "1 U | B taltn ¢ U] ‘ 7.1 fj ‘ ia“y" dy dr. (10) 


(29)? (ro+ Po)? 9 (ro+2r)(Ro—z) 


> B Ry Tro 
a ° 
2 (7° rr) Ry Ir) 
The integration over zs only extends from 0 to Ro, because the contribution from the integrand 
for x >Rp is negligible due to the presence of the rapidly varying function « 


where 


2igx 


3. Numerical Solution 


The mutual impedance Z between the two short antennas A and B for the homogeneous 


ground of electrical constants ¢ and € is given by 


h A > pw 3 Kk ’ | l 
. ” in F(ro+ Ron) | 14+. —— = * (11) 
Ir(7o ? Ro) adie u 1BU'9> Ro) B°(7'o4 Ro) 


It is now convenient to express the mutual impedance Z’ between the antennas on the inhomo- 


geneous ground as follows: 


_— h hyip 8 h ” > l l 9 
y Tt Ms. . - — —, iy (12) 
Z 27(7'o if R, } ‘ f 4 fi 7 ”) I t 1B(7'o | Ro) B(7' T E> | ; 


























where F’ is the unknown attenuation function. The latter bracketed term can be replaced by 
unity because B(ro+ Ro) is large compared to 1. Employing eq (10), (11) and (12), it follows 
that 
. . ‘ [2 "*®, F(rote2z n)F'(Ro- “£1 n)dzr 
F’(r- Ronn )~F (ro4 Ro,n) 1 Eg )( : ) ) ~ iE (13) 
“0,7 TN rip Jo (rot x)(Ro—2z)}! 
This is an integral equation for F’. It is immediately epparent when r< 0 such that A and B 
are to the left of the boundary 


F'(ro4 Ronn ~F (ro t Po,n). 


In other words, subject to our stationary phase approximation, the attenuation function F”’ 
between A and B, when they are located on one side of the boundary, is characteristic only of 
the electrical constants o and ¢ of the homogeneous ground between them. For exactly the 
same reason, the function F’(ro+2,n,m), where it occurs in the integrand of equation, for the 
general case of + >0, can be replaced by F(ro+-2,m), which is characteristic of propagation 
from A to the point r+(>0) over a homogeneous ground of electrical constants o, and «. 
After a change of variable, the final expression for /’ can be written in dimensionless form, 
vielding 
F’ =F (pp) i(”) (1—yK) pores EXP) E\pe— 39) 4, (14) 
" Jo \ PC Po Kp) 


with ° 


Fi p)=1 i2p' e” e~* da 
« ip! 
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‘a — ire re) n ) 
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“¢ > 
0; Tt tWE) , Ry 
and = | . 
a+ iwe Ro+ro 


kK 


The above expression for F’’ is then given in terms of the attenuation function F(p 9), charac- 
teristic of propagation between A and B on a homogeneous ground, and a correction term that 
accounts for the change of the electrical constants of the ground to the right of the boundary. 
It should be noted that pp is the “numerical distance” between A and B with regard to the earth 
medium to the left of the boundary, whereas poV/A is the numerical distance between the 
boundary to the point # with regard to the earth medium to the right of the boundary. Un- 
fortunately this expression, which involves an integration over products of error functions, is 
not readily expressible in closed form. However, the integral can be approximated in certain 
limiting cases. For example, when po(l1—V) and poV/K are small compared to 1, the power 
series expansion for the functions F(p) and F(py—Kp) can be employed. This has the form 


I’(p)=1—Unrp) 2p+iyrp'?s (15) 


It then follows without difficulty that 


‘ . o:{ Po ] V ‘ ; l , V v8 . 
h F(p)+2i( 2) (1-3) tan (. pr 20 \R)| & (1 J1—1 \+(z) } (16) 
For large distances, the function F must be modified to account for earth curvature Chis will be considered in section 2, 











plus terms in p,, p32, etc. Another special case is when po is large compared to 1, so that only 
the leading term in the asymptotic expansion of F'(p)— Ap) need be considered ; and furthermore, 
V is assumed to be small compared to unity. In this instance 


Fipo—Kp) _ I l 


- / 


(Po K p)' _) Po Kp) 2p 


The power series formula for /(p) is then emploved, vielding, after integrating term by term: 
i ; 7 : : i 


: + Pit - : 
) ; 4 oo 2 6 1D ¢ 


F’ ~~) —20U VK) |» _iyrp,_ 2p, ivym ,, 4D), _iyt pi__ 8 p, wr | a7 
F po) \" 3 & 7 ’ 
where p)=poV/K. 

There are probably other limiting cases which will enable the integration with respect 
to p to be effected; however, it is believed that for application to low and medium radiofre- 
quencies the values of py are in general neither large nor small compared to one. With this 
in mind, it is considered desirable to evaluate F’ by a numerical integration for a range of 
po, V, and K. It should be noted thai, for the general case, py) and A are complex so there are 
actually five parameters to consider for the two-media problem. In this paper, attention 
will be restricted to frequencies where the displacement currents in the ground can be neglected 
That is, the ratios ew/o¢ and €w/o, are assumed to be smal! compared to unity, therefore p, 

m[(Ro+7re)/A] (e9w/o) and A(=o;/c) are real quantities. This is usually justified for frequen- 
cies less than 1,000 ke for typical ground constants [10]. 

Employing the numerical values of the function F(p) and F(po— Ap). the integral in eq 
(14) is evaluated by a graphical method. The function F’(p), VA) is then plotted as a func- 
tion of pp from 0.1 to 5 for various values of A and V in figures 2 to 11. It is believed that 
the results plotted in this form can be adapted to a Jarge number of practical situations. The 
curves are not shown beyond p,=5, since this usually corresponds to higher frequencies where 
displacemeni currents are nonnegligible. 

It is of considerable interest at this stage to compare these numerical results with those 
computed using Millington’s method [3]. His empirical formula for F’’(po, V,A) in the nota- 


tion of the present paper would read 


- _ » ~ oy... [poll—V)| F(V/K) 
F'(po, VK E Pol KL) F (po) Fip(l—VyK|Ft A | (18 
Values computed from this formula for the cases A=2 and A= are indicated by © and 


@, respectively, on figures 2 to 11. The agreement is quite reasonable, and therefore further 
support is given to the validity of the Millington method, which has been known to predict 
correctly, within 1 db or so, the attenuation over certain mixed paths. It is believed, however, 
that differences are not of sufficiently small order to be neglected if fairly precise values of the 


fields in amplitude and phase are required. 
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FIGURE 12 Field strength versus distance curves for various mired paths. 
Che conductivity in millim‘os per meter to the left and right of the boundary is shown on each curve rhe experimnental results were comnunicated t 
the author privately by P. A. Field, formerly of the Canadian Broadcasting Corporation , From parametric curves , computed by Milling. 
ton method; experimental (P. A. Field curve, assuming homogeneous ground throughout (V= AK=0 
4. An Application 
To illustrate the application of these parametric curves, the field strength in millivolts per , 
meter is shown plotted in figure 12 as a function of the distance (=r,)+R)) in miles for five ¢ 
different mixed paths. The curves are normalized so that the field strength at 1 mile is 
100 mv/m. Experimental values supplied by P. A. Field are also indicated. The curves ' 
calculated by the Millington method are also shown in figure 12. The conductivities on the 
near and far sides of the boundary are indicated on the figures and are expressed in millimhos 
per meter. The agreement between the computed and experimental results is quite good. 
For these short distances the curves calculated by the Millington method fell slightly below \ 
the curves computed from the integral formula. In view of the uncertainty in the exact 
physical features of the ground, the difference between the two methods of calculation hardly 
seems significant for these situations. It is interesting, however, to note that the Millington i 
method in four out of five cases underestimates the recovery effect to a greater degree than the 
integral method. The dotted curves in figure 12 correspond to the case where the ground is 
homogeneous throughout, and has the conductivity of the earth on the transmitter side of the ; 
boundary. The difference between the solid and dotted curves is a measure of the recovery ( 
effect, and it is apparent that the theory indicates a gradual transition at the boundary and 
does not show any sudden or transient features at the boundary. It would seem from the 
experimental results that there is actually some type of disturbance at the boundary between 
the two media which is not predicted by the theory. It is believed that a more rigorous 
evaluation of the integral equation is necessary to describe the nature of the field near the 
boundary. | 


5. The Refraction Effect ” 


Although this problem was formulated for the mutual impedance between antennas A and 
B on either side of an oblique boundary, the consequent path of stationary phase did not depend 
on the inclination angle 6. In other words, any refraction effects were neglected. It is the 
purpose of this section to revise the stationary phase evaluation of the integral to account . 
for the changes of phase velocity between the two media. It is admitted that the procedure 
is not rigorous and is based on certain physical ideas that have their roots in geometrical optics. 

With regard to figure 1, the direct ray AOB between the antennas would represent the path 
of stationary phase if refraction effects were neglected. Of course, if the boundary was at a 


12 

















right angle to AOB (i. e., @=90°), it would be rigorously justified to assume AOB was the path 
of the stationary phase. It can be expected, however, that for the oblique boundary the 
path of stationary phase would be along a line APB, where P was displaced along the boundary 
from O by an amount g. Furthermore, it can be anticipated that g would be small compared 
to r, and Ry. The location of the point P can be best obtained by regarding q as a variable 
quantity in order to find when the total phase along APB is stationary or when it is a minimum. 
This is essentially a statement of Fermat’s principle in optics. The difference between the 
phase along the paths AOB and APB is now considered to be due to two factors: (1) the actual 
increase of path length, and (2) the change of the ratios of the path lengths for the two respec- 
tive media. For example, if the left-hand medium is more poorly conducting than the right- 
hand medium (i. e., A >1), the phase velocity is relatively less to the left of the boundary, so 
that AP is less than AO. More explicitly, the total phase ® over the path APB can be written 


&(q)=B(7,4 R’) +-o'(V,K, po) (19) 
where 


r rot qg°—2rog cos 6|” 


Ro=|Ri+¢?+2Roq cos 6|* 
A ‘ A a 
\ Ro ' | A BUY T Ro) Ul 
rm a 6M Po —_ 9 ( > ): 
rot Ro ™ 120% 
The function ¢(V,A,p,) in the above expression is taken to be of the same form as the function 
¢(\',K,po) previously computed. The general scheme is then to vary q to find where #(q) 
is stationary. For the present purpose, however, it is convenient to make some further 
approximations utilizing the fact that qg is small. The phase function then becomes 


| sin’ 6 /q\" l - en , 
?(q) Biro + Ro) | 14+ 9 (") va 1 OK do), (20) 
where 
P~V+! cos 0. d=rot+ Ro, 
d 
and 


“re sin’? @/q\" | ~ 
” po 1+ » (a) va 1 ni 


The terms containing higher powers in (q/d) have been neglected and, in fact, in the expression 
for py, the term containing (g/d)? can be also dropped, since ¢’(V,K,po) is a relatively slowly 


varying function of py. It is now convenient to express ¢’(1”) as a Taylor expansion, as follows: 
ii sade gq cos 6) dd’(V) qg cos 07 1] d’*o(V) 
o'(\ '( +" R A + : R 9 A (21) 
rot flo dV ro 0 — d\ 2 ~~ 


Because q is small compared to ro+ 2o, only the first term or two of the expansion are significant. 


The phase function now becomes 


hiatal q cos 6 , sin* 6 d(qy 
&(q)~B(ro + Ro) +0’ (VK, po) — 5(V)+,, ~ © ( (22 
where 6(V) do(V)/dV is essentially a positive quantity. The right-hand side of the 
equation is a quadratic in g and has a minimum or stationary value when 
qq cosé5(V) Vi—V)~A (23) 


d d . 2 sin? 6 ro Ro) 
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This value of g, namely, go, which makes #(q) stationary, is then a measure of the deviation of 
the direction of the phase of the ground wave as it crosses the boundary. The appropriate 
values of 6(V) can be obtained directly from the parametric curves of ¢’(po,V,AK). Actually. 
it is a little more meaningful to consider the angles 6 and 6,, which are the angles OBP and OAP 
in figure 1. These can be called bearing errors, as they are a measure of the change of direction 
of the phase front relative to a direct unimpeded wave between A and B. In terms of qo, the 


bearing errors, expressed in degrees, are given by 


qi _ 180 
= -- SIN 
rot+Ro O1 V) T 
and 
, Jo Pex 180 
0; > >sing 
/ to | T 


Employing typical values of the parameters and jetting A= © , values of 6 and 6; were com- 


puted and are: 


@—445 6é—60 
| Ey R 
o \ 
] c Oo i) 

0. 2 0. O1 50 0. OG 0. 24 0. 03 0. 14 
2 Ol 100 03 14 . 02 . OS 

2 OO] 50 02 O08 01 O4 

2 OO! 100 Ol O6 01 03 
5 01 50 _10 . 10 . 06 . 06 

5 01 100 . 06 . 06 . 03 08 

5 OO] 50 08 03 0? 02 

5 OO! 100 02 . 02 Ol O01 

8 Ol 50 26 06 .17 04 
8 01 L100 . 16 04 . 09 02 
= OO! 50 OS . 02 O5 Ol 

8 001 100 . 06 . 02 . 04 . 01 


It might seem surprising at first glance that these values are so small. This is as should 
be expected, however, because the absolute phase velocities of the ground wave over the two 
media differ only by a fraction of one percent [10j. It is noted that the bearing errors are 
largest when the antenna is near the boundary, and where the separation between antennas 
A and B is 50 wavelengths or less. Of course, as the obliqueness of the boundary becomes 
more noticeable (i. e., 045°) the calculated values of 6 and 6, would be larger. One should 
be cautious, however, in attempting to apply this method to the case where the boundary 
makes a small angle with the line AB between the antennas. In such an instance the approxi- 
mation used in evaluating the basic integral by a stationary phase principle becomes invalid. 


6. Conclusion 
The results obtained in this paper provide a theoretical check of the Millington semi- 
empirical method of calculating the amplitude and phase of a wave crossing a boundary sepa- 
rating two media. The parametric curves presented here should be convenient for making 
predictions for propagation along a two-media path at low and medium frequencies. It is 
also indicated that the refraction effects for a wave crossing the boundary obliquely are very 


small. 


I record my thanks to Loris Perry, who carried out most of the computations. 
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Effect of Object Frequency on Focal Position of Four 
Photographic Objectives 


Fred W. Rosberry 


Three wide-angle aerial-camera lenses of 6-inch nominal focal length and one smaller 


lens of 32-millimeter focal length were given photographic resolving-power tests. 


Targets 


were of the parallel line type with several groups of three lines representing a wide range of 


lines-per-millimeter values. 


The image plane of best focus was visually located, and expo- 


sures were made at four different positions inside of this focus and the same number outside. 
Measurement of the resultant negatives by a recording microdensitometer revealed the image 


contrast of the various line-space frequency groups at the nine focal positions. 


Curves are 


presented showing the focal positions of maximum image contrast for four values of line-space 


frequency. 


1. Introduction 


During the past decade, there has been consider- 
able interest in the effect of target contrast on image 
contrast and resolving power in the image plane. 
Investigation of the relationships between these 
quantities has shown the existence of a number of 
associated phenomena, among which is a shift in 
the plane of best focus with target pattern frequency. 
Any appreciable shift will affect the selection of the 
proper focal plane in a fixed-focus camera, such as 
an aerial camera [1, 2].'. The Optical Instruments 
laboratory has therefore made an investigation of 
this phenomenon. 

Of the lenses on hand at the time this investigation 
was being planned, three, representing a good, 
medium, and fair quality as shown by previous 
resolving power tests, were chosen for testing. They 
were all wide-angle aerial camera lenses of nominally 
é-in. focal length, made by the same manufacturer. 
A fourth lens was later included representing a 
diferent tvpe and of only 32-mm (1.3-in.) focal 
length. 


2. Apparatus and Procedure 


To measure the contrast rendition, or contrast 
in the image plane of a lens, by a photographic 
procedure, extreme care had to be taken to control 
the other factors that can affect contrast. Some of 
these factors are exposure, development, and varia- 
tion in emulsion from one photographic plate to 
amother. Uniform exposure was maintained by the 
ue Of an automatic timing device, with relatively 
long exposure times such as 2 or 3 sec instead of 
small fractions of a second as obtained by shutters. 
Development was maintained uniform by putting as 
many different exposures as possible on one plate and 
by agitation during development. The variation in 
characteristics between plates out of the same box 
Svery small, but this was further reduced by cutting 
large plates into smaller ones, thus making two or 


more from one. 
eeeeeeeeeeenns 
Figures in brackets indicate the literature references at the end of this paper 
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The targets used in this investigation were the 
75-line step contrast and the 200-line step contrast 
targets, which are variations of those previously 
reported [3]. The 75-line target consists of 25 groups 
of 3-line patterns. The line widths and spacings in 
a given pattern are equal; however, the values vary 
from one pattern to the next by a ratio equal to 12. 
The renge of resolving power for this master target is 
therefore from 0.2 to 3.2 lines/mm. The range of 
the 200-line target is from 0.2 to 2 lines/mm, which 
means that the actual widths of lines and spaces vary 
from 2.5 to 0.25 mm. This range is covered in a 
continuously varying linear manner with 200 lines 
and 200 spaces contained in a distance of 180 mm. 
The zero line is 20 mm from the origin, so that the 
over-all length of the resolving-power portion of the 
chart is 200 mm (fig. 1). The length of the lines is 
185 mm, thus giving the contrast scale approxi- 
mately the same length as the frequency scale. 
Along the length of the lines the targets are divided 
into nine areas of different line-space contrast values. 
The 75-line target had a contrast range in these nine 
areas from 0.22 to a high of 3.19 on the density scale. 
The same nine areas on the 200-line target had a 
range from 0.10 to 1.51. 

Although two targets were used in this investiga- 
tion, they were not used simultaneously. The target 
under observation was mounted on the axis and in 
front of the lens being tested. As the targets are 
transparencies, it was necessary for them to be uni- 
formly illuminated from behind. This was accom- 
plished by using a light box of such design that only 
diffusely reflected light, which was quite uniform and 
free of hot spots, fell upon the target. 

The camera used for exposing the test plates 
consists of two major components, namely, the lens 
mount and the camera back (fig. 2). A set of ways 
similar to a lathe bed supports both components, 
each mounted on its own carrier or saddle, thus 
allowing for the necessary movement of one with 
respect to the other by means of a lead screw. 
This movement is necessary to accommodate lenses 
of different focal lengths. The camera back itself 
is provided with three micrometer screws, by means 
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Figure 1. The 200-line step contrast target as used in this investigation 
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Figure 2. Bench showing lens, plate carrier, and focusing 
micros ope 
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af which small precisely measured axial displace- 
ments of a millimeter or less may be made. The 
plate holder can also be moved vertically In & plane 
perpendicular to the lens axis, allowing successive 
images to be taken on a single plate, one above the 
other. Changes 1n the foeal position of amounts as 
small as 0.2 mm are easily accomplished with the 
micrometer adjustment. A microscope is mounted 
on the same saddle as the camera back and focused 
on the plane of the emulsion. 

The lens mount supports the lens rigidly in place 
and carries the necessary bellows and light shields 
to protect the negative from stray light. The mount 
constructed so that it holds a metal lens board 
measuring 6 in. square, which in turn carries the lens 
under test. 

The distance between lens and target was equal to 
% focal lengths of the Jens under test. With this 
lens-to-target spacing the image distance is a little 
longer than for an iniinitely distant target, but it is 
believed that there is very little change in the resolv- 
ing power characteristics. Theoretical considerations 
indicate that the lens-to-target distance can be re- 
duced to as low as 15 focal lengths without signifi- 
eant change in the resolution characteristics from 
those for infinitely distant objects.’ 

With the lens-target distance set at 26 focal lengths 
for the 6-in. lenses, the effective spacing of lines in 
the target on the image plane was 5.0 to 79.6 lines 
mm for the 75-line target. With the 32-mm lens, 
the lens-to-target distance was 51 focal lengths, 
giving a range of 10 to 100 lines/mm for the 200-line 
target. Focusing was accomplished by properly 
spacing the camera-back and the lens until a well 
defined image appeared on a ground glass, which was 
in the position the photographic plate would occupy 
during photographing. The microscope was then 
focused on the plane of the ground glass and locked 
in that position. With the ground glass removed, 
it was then possible to view the image directly and 
make any small focusing adjustments necessary. 
The micrometer was then set so that the plate carrier 
was 0.8 mm inside of the best visual focus, and a 
panchromatic photographic plate set in the carrier. 
Using the vertical ways of the camera back, the plate 
was racked to its uppermost position and the first 
exposure made. The plate was then moved away 
from the lens by 0.2 mm and lowered enough to 
allow the image to be put above the first and another 
exposure made. This procedure was repeated until 
nine exposures were made, covering the range from 
0.8 mm inside to 0.8 mm outside of the visually 
estimated focus. After developing in Dk-60a for 5 
min and fixing, this plate was then ready for measure- 
ment. 

The density of each line and space for three of the 
target contrast levels were the quantities to be 
measured. These were determined by using a re- 
cording microdensitometer. The test negative was 


—$________. 


? Private communication from D. P. Feder 
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located on the plate carrier of the instrument, prop- 
erly oriented, and brought into good focus. One 
contrast area of one focal position of the test negative 
was then scanned and the local values of line and 
space density were automatically recorded in the 
form of a graph. The negative was then relocated 
so that the same contrast area on the next adjacent 
image, representing the next focal position, was then 
ready for scanning. The instrument was refocused 
at each change in negative position. All focal 
positions were covered by this procedure, and this 
constituted a run. Other runs were made for the 
two remaining contrast areas to be measured. Each 
run was repeated two more times with the results 
averaged to minimize any variation in measurement 
that still might exist, thus making three runs of the 
three contrast areas, with each run covering the 
nine focal positions. 


3. Results of Measurement 


The graph produced by the microdensitometer con- 
sisted of curves on a strip of chart paper graduated 
with a logarithmic scale of photographic density 
whose limits were zero and infinity. The quantity 
measured was the amount of light transmitted by a 
very small area of the test negative. Because of the 
line and space characteristic of the negative being 
scanned, the recording pen was traveling from near 
one end of the scale to near the other, thus indicating 
the highest transmission or minimum density at a 
space and the maximum density at a line. These 
data were tabulated by reading, on the log scale, the 
naximum and minimum tor each cycle. The dif- 
ference of the two was called contrast (C). The 
75-line target, which was divided into 25 groups of 3 
lines and had the line frequency constant for each 
group, was averaged to determine the maximums and 
minimums for each group. A graphical representa- 
tion of these data is shown in figures 3 and 4. The 
curves are grouped according to the particular lens 
under investigation and the target contrast. Each 
curve represents the amount of contrast in the image 
at various focal positions for a specified image fre- 
quency. Some of the curves show a bar indicating 
the calculated depth of focus for the adjacent fre- 
quency or resolving power. The depths of focus 
herein referred to were calculated from the formula 
d=4b/a, where d is the depth of focus in millimeters, 
b is the f-number, and a is the resolving power in 
lines per millimeter, for which the depth of focus is 
sought [3]. This formula was developed for use with 
high-contrast targets, which should be taken into 
account when viewing the results. Theoretically the 
curve should reach zero value of image contrast at 
about the same points as indicated by the ends of the 
depth-of-focus line. The range of the depth of focus 
is shown only for the two curves of highest frequency, 
as the range of focal positions is not great enough to 
include the depth-of-focus range at the lower 
frequencies. 
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Fictre 3 Curves showing the contrast in the image at various focal positions of fow line-fre quency patterns at three values of 
target contrast for the three wide-angle aerial camera lenses. 
L/mm, lines per millimeter; Di, photographic density of line; ),, photographic density of space, Che lines marked “DD” indicate the range of the calculated dept 
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4. Discussion 


Contrast can be expressed as the difference be- 
tween two densities (A//), or it can be expressed as a 
atio of the two transmissions T/T. Both methods 
ve in use, with perhaps the former more common. 
some of the curves were plotted by both methods, 
and although each illustrated the point well, the AD 
procedure lent itself better to data from the micro- 
densitometer. One point rather clearly shown by 
se curves the relatively small amount of 
suift in the optimum focus for the four frequencies 
blotted. This of course is a quality that will change 
irom one lens to another and is somewhat dependent 
mi the individual lens aberrations 

The three wide-angle aerial camera lenses were 
given resolution tests prior to this investigation, and 
the results of these tests established that their reso- 
ution characteristics were within the range covered 
by the average lens of this ty pe. Distortion measure- 
ents also indicate that the sample lenses are not 
outside the average. The lenses tested herein are 
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Fieure 5. Curves showing the shift from an arbitrary zero of 
the focal position giving highest contrast for various line 


frequencies at two target-contrast levels. 


The same three lenses are shown in figure 4. 


quite consistent in showing very little shift of focus 
at the two higher target contrasts; at the low target 
or object contrast (C,—0.22) the shift is even less 
noticeable. The sets of curves shown in figure 5 
summarize the findings for the three 6-in. wide- 
angle lenses. The amount of shift (from an arbi- 
trary zero) of the focal position giving highest con- 
trast is plotted for four values of image frequency. 
The two plots represent the results obtained with 
two different target contrasts. The focal position 
measured along the ordinate of the graph was de- 
termined from the crest of the corresponding con- 
trast curve shown in figure 3. The contrast curves 
for the lower frequencies are rather flat-topped, 
which makes precise location of the peaks difficult. 
Shifting of the plane of best focus with object fre- 
quency as indicated herein is very small for lenses 
of near average quality. 

As focal settings become more precise, the position 
of best focus becomes more difficult to locate. <A 
zone is found within which a position of best resolu- 
tion and another position of highest contrast rendi- 
tion are located. This effect is present at all points 
in the field and, as the range of this zone of best focus 
is less than the inherent field curvature of a lens, 
critical focusing over the entire field becomes ex- 
tremely difficult. The residual spherical aberration 
in the lens is probably a cause of this focal shift. 
The position for maximum contrast of broad lines is 
where the circle of confusion is the smallest and not 
necessarily where the resolution is best. 

The present work accordingly confirms the results 
of earlier workers [1, 2, 5] that a difference in focal 
setting exists between the settings for best resolution 
and for best contrast rendition. However, it em- 
phasizes the fact that this difference is usually so 
small that it is seldom advantageous to set for the 
position of maximum contrast rendition for low- 
frequency patterns as opposed to setting for the 
position of maximum resolution. For example, in 
some earlier work [1, 2], the difference in focal setting 








for frequencies in the image of 4.1 and 20.6 lines/mm 
was indicated to be as high as 0.03 in., or 0.76 mm. 
It was further indicated that the improvement in 
contrast rendition would be quite appreciabie at the 
lower frequencies in setting at the optimum image 
contrast for those frequencies with marked loss at 
the higher frequencies. It is clear from the results 
shown in this paper for the three lenses, shown in 
figure 3, that the shift in focal settings between 
frequencies of 5 and 20 lines/mm is at most 0.20 mm. 
Moreover, the slight gain in contrast rendition by 
setting at the maximum for the low frequency is not 
sufficient to warrant accepting the appreciable loss 
in contrast rendition at the higher frequencies. It 
is probable that some few lenses of a given type show 
focal shifts as bigh as 0.76 mm, but it seems likely 
that such lenses are the exception rather than the 
rule 
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An Examination of the 1955 Helium Vapor-Pressure 
Scales of Temperature 





E. Ambler and R. P. Hudson 


In a previous communication, magnetic thermometer calibrations in the region 1.3° to 
1.2° I were used to examine the internal consistency of two provisional helium vapor-pressure 
scales of temperature. Since that article went to press, these scales have been made avail- 
able in their final form, and further measurements have also been made. The results of all 
the measurements, seven runs in all and with two different apparatuses, are used here to 
examine the latter scales. The results of six runs were self-consistent and could be reconciled 
with the empirical scale of Clement et al. (7 ;—.) towithin 2 millidegrees, and with the cal- 
culated seale of Van Dijk and Durieux (7;;) to within 5 millidegrees. The remaining run 
showed better agreement with 7;;. It appears, therefore, that the previous assessment of 
the 7; scale may have been a little too favorable. The most critical interpretation of the 
findings is that the present technique of vapor-pressure measurement is not sufficiently 
precise to discriminate between the two scales, 


1. Introduction in very close accord with those obtained with the 

“first apparatus.” 
In a recent publication [1]! the authors examined The present paper deals with an examination for 
the internal consistency of two vapor pressure- | internal consistency of the 7;; and 755, tabulation, 


ibsolute temperature (p-7) relations through the | based upon all of the magnetic thermometer calibra- 
calibration of a magnetic thermometer in the range | tions referred to above. 





13° to 4.2° K. One such relation was that of 


| Clement, Logan, and Gaffney [2], based upon an 2. Experimental Procedure 

empirical equation, and the other a thermodynamic 

alculation by Van Dijk and Durieux [3}. | The apparatus has been described in detail in [1] 
After reference [1] had been prepared for publica- | (see fig. 1, b, and text therein). The only modifica- 





tion both the new p-7' relations were modified by | tion in procedure was the employment of a mano- 
their respective authors. Clement based his final | stat ? to control the temperature of the bath. This 
temperature scale [4] (to be referred to by 75x) on | proved very useful in that the system came to equilib- 
the experimental data of several authors, including | rium twice as fast as previously and enabled one to 
imbler and Hudson, obtained in one of two ways, | obtain twice as many points per run. (It may be of 
iz, a measurement of the vapor pressure using (a) & | jnterest to record that at the lowest temperature 
separate vapor-pressure bulb, or (b) the pressure employed (1.3° K), where the magnetic thermometer 
wer the liquid helium bath plus an appropriate | js extremely sensitive, the bridge reading remained 
orection for the “hydrostatic head,” when a heat | constant to within the equivalent of 107° deg for 
source Was employed to maintain a constant flow of periods in excess of 5 min.) 

tubbles up through the liquid. Van Dijk [5] has 

nereased his temperature values slightly in the region 

3° to 4.3° K . that his p-7' curve is tangent to a 3. Results 
rve through the gas-thermometer data of Berman 
ind Swenson [6] near the normal boiling point. 
Temperatures derived from this new seale will be 


The data obtained in the new series of measure- 
ments are summarized in table 1. The pressures, 
—~ p, have been corrected to 20° C and standard gravity, 
lenoted by l : and the bridge readings, n, have been corrected for 
Concurrent with this work, the present authors nonlinearicy in the mutual-inductance decades. 
arried out further measurements, employing the The data of table 1 and those of table 1 in [1] will 
second apparatus” of [1] to augment the meager | now be discussed in relation to 7;; and Tsp. 
mount of data obtained therewith. This was 
specially important as these same data were not 


The authors are indebted to L. D. Robertsofthe Oak Ridge National Labora- 
— tory for a drawing of this device. The original design is due to H. 8. Sommers, 
figures in brackets indicate the literature references at the end of this paper. Los Alamos, Rev. Sci. Instr. 25, 793 (1954 
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TAPLE 1. Corre sponding values of pressure in millimeters of mercury, p (corrected to 20° C and standard qravity) and bridge dig 
setting, n (corrected for nonlinearity in decade scale ? 
January 4, i956 January 18, 195¢ 
p ’ D p n p p i] 
751. 42 17. O77 124. 33 20. 5650 0. 46 17. 7534 4). 395 22. 854 13. 398 25. 1098 
717.24 17. 6786 17. 26 21.0243 670. 452 17. S068 7. 766 23. 0118 9. 016 25.9244 
716. 68 17. 6797 7. 004 22. S817 583. 62 18. 1039 $2. 392 25. 3254 5. 968 26. 7415 
628. 73 17. S697 1.771 23. 2392 184. OS 18. 3710 32. 418 23. 3208 4.017 27. 6820 
506. 19 18. 1953 24. 846 23. 7329 0). 65 IS. 8405 24. O12 3. W254 1.031 27. 6271 
357. 75 18. 7391 14.045 24. 7613 197. 36 14. S600 24.115 235. % 168 >. US4 28. 2729 
241.0 19. 3912 &. 123 26. O191 101. 97 21. 0572 18. 124 24. 49358 633 24. 6006 
3.1. Comparison With Tss Temperature Scale TaBie 3. Values of constants A and B in equation 
i | evaluation 
Che data were analyzed according to the procedure 
outlined in [1]. When the constants A and B in Date K 
the formula 
June 7, 195 34. 1649 19. 3198 
, al June 8, 1955 34. 1461 19. 3027 
n BR A ] | June 16, 1955 $4 1703 19. 3187 
June 17, 1955 $4. 1555 19.310 
(August 17, 195 23. (i389 11. 327. 
have been determined, a value of temperature, ey a ee ie ae po 
— nual , th 23. H962 2 09 
“magnetic temperature,” 7, is at once calculable 
for a given value of n. The corresponding value of p 
leads to 7;, from the p-7 table and AT=T. T, ‘ , 
I i comparing with table 2 With this temperatw 


mav be plotte “las a function of absolute temperature 
This is done in figure 1. The values of A and B used 
for each run are summarized in table 2 


rarLeE 2 Values of constants A and Bin equation 1 for T 
eval sation 
Date I 
June 7, 1055 4 O56 G 284 
June &, 105 4.008 14. 2644 
June 16, 19 +0 14. Jue 
June 17, 1955 4 O44 19, 2809 
August 17, 1955 93. 5T6 l $537 
January 4, 1956 23. HO73 W922 
January 18, 105¢ 23 OHO 9 115 


is evident that only the results of August 17 


lt 


1955—the first run with the second apparatus —can 
be brought into very close accord with the 7 
tabulation The rest of the data can, however, be 
brought into very close harmony with each other 
by suitable choices for A and PB in each case, and 
the entire group deviates considerably from the 
zero line with a maximum of +4.5 millidegrees at 
3.6° K. There appears, furthermore, to be a dis- 
continuity in the region of the lambda-point. (As 
deseribed detail in [1], the choice of A and P is 
made so as to give an optimum fit throughout the 


The implicit assump- 
by a 


entire range of measurement. 
tion is that the 7. scale is nowhere in error 
large amount. 


3.2. Comparison With Tsszs Temperature Scale 


Figure 2 was obtained by a procedure analogous 


to that leading to figure 1, now using 7;,;. instead 
of 7;;. Small changes in the values of A and B were 


hecessary as May be seen by consulting table 3 and 


scale it is possible to restrict the value of AT with 
the limits of +2 millidegrees. The discontinui 
at the lambda-point is still evident. (Although su 
a discontinuity is more likely due to faults in tj 
measuring technique than genuinely anomaloy 
behavior of the p-7T curve, it has been observed | 
other investigators. See, for example, Erickson an 


Roberts [7]. 


4. Discussion 


ln [1] reasons were given for 
in the results obtained with the second apparatw 
(August 1955 data) than in those for the first appr 
ratus (June 1955). These were based upon the a 
tempted improvement in design and the fact that: 
rather rough) measurement of the bath pressun 
plus the hydrostatic-head correction agreed with | 


measured in the bulb svstem, within the limits ¢ 
error of measurement. It can be seen, however 
that the measurements of January 1956 reprodw 


the earher results obtained with a different apparatus 
and fail to agree with the first measurements, 
the same apparatus. 


mentioned criteria, had not changed. We have be 
unable to find a plausible explanation for the dis 





using 
At the same time, the quality] 
of the second apparatus, as determined by the above 


reposing more trust} 


‘| 


| 
| 


crepancies, which must therefore be held to indicat 
the degree of trustworthiness for the entire series ¢ 


experiments. Hence, while the present results viewe 
as a whole favor 7.» to the disadvantage of 7 
it is probably correct to conclude that the present 


used techniques of vapor-pressure 


stant .1 decreased slig 


ery small, seems! 


the con 


rT Tect. 1 


It may be noticed, further, 
on the second day of a 2+lay serie hough 
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ge dig} yot sufficiently precise to resolve the discrepancy 5. References and Notes 
between the two scales. This latter is illustrated in : 
: a Rais Smee lotted T T.. : ; [1] E. Ambler and R. P. Hudson, J. Research NBS 56, 99 
figure 2 where we rave plotter 55 5,5 Us a (1956) RP2654. 
function of T. [2] J. R. Clement, J. K. Logan, and J. Gaffney, Phys. Rev. 
‘ } 100, 743 (1955). 
In conclusion, it should perbaps be emphasized (3) H. Van Dijk and M. Durieux, Proe. Conf. de Phys. 
. . . . P — — 28 Pa rio Pre. y 
that although use was made of the results in [1] in oa ~~ remperatures, Paris 1955; Commun. No. 
(aes , oe 7 0, Pp. oo. 
the development of Tss5, no greater weight was | [4] J. R. Clement (unpublished). This seale is the one 
seorded them than any other of four additional referred to in the “note added in proof” to reference 2. 
e ; - | ; [5] H. Van Dijk (private communication). It is understood 
> » ré =~ Sh " ) ’ " ys) ; . . 
independent investigations over the same temp¢ von that copies of this scale are available on request to the 
ture range. This accounts for the systematic devia- Kamerlingh Onnes Laboratory. 
Vor olla ” 9 9°1g 9 gO JQ0T [6] R. Berman and C. A. Swenson, Phys. Rev. 95, 311 (1954). 
» Tre Ss : Z wnd o5.¢ O 4.2 7 "_e : = 
tion in the regions L.3° to 2. K and 3.6° t . K 7] R. A. Eriekson and L. D. Roberts, Phys. Rev. 93, 957 
visible in figure 2. (1954) 
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Heats of Formation of Xonotlite, Hillebrandite, and 
Foshagite 


Edwin S. Newman 


The heats of solution of synthetic xonotlite, 5CaO-58i0O,-H,O, hillebrandite, 2CaQ-SiO,- 
H.O, and a sample of naturally occurring foshagite, 5CaO-3Si0O,.-3H,O, were determined in a 
mixture of nitrie and hydrofluoric acids, and their heats of formation from the oxides were 


92.6, 34.7, and 


calculated to be 


1. Introduction 


The system CaQ-SiO,-Al,O,-Fe,O;-H,O includes 
the major constituents of portland and high-alumina 
cements and the reaction products formed by the 
hydration of these materials. In connection with 
the hydrothermal study of parts of that system 
1.2)! being conducted at the National Bureau of 
Standards, the heats of solution of the reaction 
products and of related minerals are being deter- 
mined. The heats of solution of these materials and 
of the reactants from which they theoretically can 
be formed are measured in appropriate acid solutions 
in a precision-type calorimeter. The data, com- 
bined with literature values for the heats of forma- 
tion of the reactants [3], permit the calculation of 
the heats of formation of substances related to the 
products of hydration of the cements. Appropriate 
thermochemical equations can be written and the 
computed heats of reaction can be compared with 
observed measurements [4,5] for the further elucida- 
tion of the hydration reactions occurring between 
water and hydraulic cements. The present paper 
gives the results of measurements made with calcium 
hydroxide and silica gel as reactants and of synthetic 
hillebrandite, 2CaO-SiO,-H,O, and xonotlite, 5CaQ- 
js10,-H,O, as products. In addition, the heat of 
lution of a sample of the mineral foshagite was 
determined. 


2. Materials, Apparatus, and Procedure 
2.1. Materials 
Reagent-quality nitric and hydrofluoric acids, 
calcium carbonates (low alkali), calcium hydroxide, 
and precipitated silica (silicic acid) were used. Cal- 
cum hydroxide for heat-of-solution measurements 
was prepared by igniting CaCO, for 72 hours at 
1,000° C and treating the resulting CaO with water 
10} at 250° C for 2 days. The crystals of Ca(OH), 
obtained were dried, without washing, in a vacuum 
over magnesium perchlorate and stored in a tightly 
stoppered bottle in a desiccator over saturated KOH 
mtil used. The ignition loss of this material was 


rackets indicate the literature references at the end of this paper 
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94.6 kilocalories per mole, respectively. 


24.30; percent, to be compared with 24.31, percent 
calculated from the formula. 

The reagent-quality silicic acid used for the heat- 
of-solution measurements had a water content of 
12.82 percent and a nonvolatile residue of 0.023 
percent when treated with hydrofluoric acid. The 
silica was kept over saturated MgCl,-6H,O solution 


for 2 weeks before its heat of solution was deter 
mined. At that time the water content was about 
20 percent. According to Mulert [9], whose data 


were used in calculating the heat of formation of 
SiQ,-aq given in reference [3], the heat of solution 
of hydrated silicic acid per mole of silica does not 
change with further increase in water content. 

Xonotlite was prepared from Ca(OH), and silicic 
acid. These materials were mixed together in a 
CaO/SiO, ratio of 1.000 by tumbling for several 
hours in a glass bottle with a few wooden balls. 
Twenty grams of this mixture was prepared and 
heated overnight at 950° C in a platinum dish, 
cooled in a desiccator and flooded with freshly boiled 
distilled water. The dish was transferred to a small 
bomb and heated at 250° C. The flooding and 
transfer of the dish and the closure of the bomb were 
made as rapidly as possible to minimize exposure to 
CQO, in the air. The bomb was opened at, the end of 
2 weeks and the hardened mass of surface-dry 
material ground quickly to pass through a No. 28 
sieve and dried for 2 hours in a vacuum over Mg 
(C1O,)>o. 

Three samples of synthetic hillebrandite, C,SH(B) 
[11], were prepared in a similar manner by heating 
the material with water at 242°, 178°, and 152° C for 
7, 16, and 14 days, respectively. These samples 
also, as well as selected fibers of foshagite from Crest- 
more, Calif., were crushed rapidly to pass a No. 28 
sieve and dried. 

The samples were analyzed for SiQ,, Fe,O,+- Al,Os, 
CaQ, and MgO by the methods described in Federal 


Specification SS-C-158e. Loss on ignition was 
determined by heating 1-g samples at 1,200° C. 


The CO, content was determined on 2-g samples, 
and the water content was taken as the difference 
between the loss on ignition and the CO, content. 
X-ray patterns, obtained by using a recording 
X-ray diffractometer with copper radiation 1.5418 A, 
were used to confirm the identity of the samples by 
comparison with published patterns [12]. 








TAPLE 1 tnalyses of calorimete) sam ple g 
XNonotlit Hillebrandite Foshagité 
Determinations Ca(OH SiOy-aq Observed 
Com : Com Com 
puted Observed | puted puted Observe 
l 2 $ i 
Ignition residue 75. 69 87.18 oH. 99 6. 54 OO. 53 67. 57 SY. 25 SY. 20 SU. 48 91.74 gO. 50 RS 4 
CaO 75. 69 45.83 45. 91 Sa. U6 13 7. 52 7. OD os. 4 60.03 4. 4 418. 41 
SiO 87.16 16 18. 8 $1. 57 23. 80 31.56 31.55 41.02 41. 86 $5.01 36. 6 
CoO os 0. 27 1. 27 1. 6S 1. 6S 0. 35 0.04 7 
H.O «4 24.16 12. 82 $01 $19 47 31. 16 ¥. 07 9. 03 10.17 8. 22 10. *) WRN 
R,O 0.21 0.09 0.12 0.12 “14 0. 26 ON 
MeO $1 4 ON ) 
HF residue 0. 023 
Analvsis, total 100. 00 100.00 100. 00 US. 7 loo. oOo un mY. Go. Y. OL ol 100. 49 100, 00 100 17 
As taken from stock botth Computed from formula 5CaQ-58iO,-H,O for xonotlite: 2CaO-SiO»,-H.O for hillebrandit« ind 5CaO-3S8iO»3HLO for foshagit 
By difference 1 Loss on ignition minus CO 
lhe results of the chemical analyses are given in | heat of solution of silica in hydrofluoric acid becomes 
table 1. For easy comparison there are also included | constant [9]. It was observed that 3-g samples of 


the analyses of the silicates computed from the 
formulas. The first sample of hillebrandite, pre- 
pared considerably before the others, was found to 
contain excess moisture. It was dried further and 
analyzed and its heat of solution again measured. 
Subsequently, its ignition loss and heat of solution 
were redetermined. This material appears as hille- 
brandite samples 1, 2, and 3 in the tables. Figures 
in parentheses in table 1 represent its analysis as 
sample 2 computed to its moisture content as samples 
l and 3 

The agreement between observed and calculated 
compositions was satisfactory for the synthetic 
preparations. The foshagite, however, had under- 
gone considerable alteration. The sample consisted 
of selected fibers, taken from a 50-g fibrous mass 
Perhaps 30 percent of the sample was of somewhat 
lower index than reported for foshagite, and 
brownish material was present. 


some 


2.2. Apparatus and Procedure 


The calorimeter, has been adequately described 
elsewhere I6, 7}. It was operated so that the final 
temperature was always within a few hundredths of 
a degree of 25°C. Since the heat capacities of the 
reactants were used in computations, the isothermal! 
heats of solution at 25°C were obtained. The energy 
equivalent of the calorimeter was determined elec- 


trically, using the defined calorie of 4.1840 absolute 
joules. The heat capacities of the samples, intro- 
duced at room temperature, were taken as 0.2 


cal/deg-g. The calculated correction for the sensible 
heat introduced with a sample rarely amounted to as 
much as 0.05 cal/g. 


The acid charge used had a total weight of 600.0 g, 
of which 11.0 ml (12.6 g) was 48 percent HF, the 
remainder being 2.000 N HNO, (at 25° C). Its com- 
position therefore IS! was O.30HF, 1.10 HNO,, 


29.10H,0. The sample weight was generally 1 g. 
When silica gel containing 20 percent of HO was 
being dissolved, the HF/SiO, molar ratio was 22.8, 


well in excess of the value of six, beyond which the 
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portland cement containing about 1.8 ¢ of combined 
CaO could be dissolved without precipitation of 
CaF,, but that when dissolving Ca(OH), alone thy 
sample weight had to be reduced to 0.7 g. Unde 
these conditions there was no weighable or visible } 


precipitate of CaF,, but with 0.8-g samples of 
Ca(OH), a faintly visible haze remained in tly 
solution after filtration had removed 0.6 mg of 


CaF,. Consequently 
0.7-g samples of Ca(OH), were used to determin 
the heat of solution of that substance. It was 
assumed that the heat of the dilution effects caused 


precipitated approximately 


by the differences in the concentrations of the fing 
solutions were negligible, and no corrections wer 
made. 





3. Results and Discussion 


The observed heats of solution are given in table? | 
The amount of CaCO, equivalent to the CQ, was 
calculated for each sample, and its heat of solutior 
86 cal g [4], deducted from the observed heat af 
solution. The total CaO was reduced by an equiy- 
alent amount. The heats of solution of the samples 
per mole of SiQ,, A, were calculated by the formula 
A a(100 60.06/P). where a is the observed heat of 
solution of the sample, corrected for CQ,, P is thi 
percentage of silica in the sample as given in table 1, 
and 60.06 is taken as the molecular weight of SiO; 
These calculated values are given in table 3. Enm- 
pirical formulas, calculated from table 1 corrected 
for CO,, are also given in table 3. In making th 
foregoing calculations for the foshagite, an equiv- 
alent amount of CaO was substituted for the MgO 
found by analysis. 

The heat of solution of C/S? moles of Ca(OH 
calculated for each sample and added to the heat ol 
solution of 1 mole of SiO,, 0.828H,O to determine the 
solution of the reactants from which the 


| 
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SH 


10, SiOe, and HeO are represented by the let 


+ SiO>oHeO heir represented by the formula ¢ 
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TAPLE 2. 


Calorimeter 
energy 
equivalen Ca(OH SiO.aq Xonotlite 
caliohm 4 cal'g caliq cal/q 
5723. 2 $29. 97 468.72 490.8 
572. 2 $24. HY 468. 64 4. 92 
5716.4 $20. 32 166. 51 490. 67 
S718. 8 429. 31 470. 86 491. 66 
5717.0 $20. 26 470.00 $00.16 
430. 24 468. 51 400), 42 
Average 5719.2 $24.63 168. G2 100.78 
Std. dev 
of ave 1.2 07 0 7 0. 21 
y sample in 600-2 HF+HN Os; mixture (0.30HF, L1IOHN Os, 29.10H2O 


0.7-¢ sample in 600-g HF+HN Oy; mixture (0.30HF, 1.10HNOs, 

Sample stored over saturated solution of MgCh-6H2O 
310», 0.02 percent of HF residue, 19.89 percent of H,O. H/S=0.82s8 

ARIAT is 0.1010 at 25° C; 5719.2 cal/ohm is approximately 577.6 cal/deg 

1-g sample in 425-g HF+HN Oy, mixture (0.22HF, 0.78HN Os, 20.61H2,O 


29.10H,O 


Results of calorimetric measurements 


Heats * of solution 


Hillebrandite 


Natural 
foshagite 
] 2 3 4 

cal qd cal g cal g cal gd cal g cal g 
* 354.4 462. 40 1). 62 470. 47 490. 39 437.94 
351.0 461. 57 459. 00 470. 32 189. 19 438. 67 
461.79 459. 53 469. 04 488. 47 438. 66 

459. 46 469.14 488. 51 439. 31 

459. 96 468. 91 488. 48 438. 67 

159. 36 469. 29 489. 39 437.27 

352. 7 51. 92 459. 66 469. 68 489. 07 438. 42 
1.7 0. 25 0. 23 0. 27 0.31 0. 29 


Relative humidity 33 percent 


Ignition residue of sample 80.11 percent, equivalent to 80.09 percent of 


in vacuum-flask calorimeter 


Standard deviation of average y{z z nijn(n—l 
TAPLE 3. Heat of reaction of Ca(OH), and SiO, O.83H.O 
Heat of solution 
Heat of 
Sample Empirical Formula reaction, 
Products Reactants AH 
keal/ mole keal/ mole keal/mole 
( OH 1.00Ca0-1.00HLO 31.88 
SiO»aq 1.008i02-0.83H2O 35. 16 
Xonotlit 0.98C aO0-1.00S8i10e-0.21H2O 50. 06 6. 40 7.34 
Hillebrand 
1 L.88CaO-1.0OSIO»-4.35H,0 8S. 4 95. 09 7.05 
? L.&&CaO-1L0OS8SI0.-0.96H2O 87. 31 95. 09 7.78 
L.&&CaO-1.00S8i0.-0.96H2O 86, 82 95. 09 &, 27 
ry 2 01Ca0-1.008i02-1.09H,O oO. 80 OY. 24 8. 44 
2 02Ca0-1.008i100-0. 86H2O 92.18 09. 56 7. 38 
Weighted average 7.91 
Foshag 1.44Ca0.-1,008i102-0.81H2O0 4 71.13 81.06 9.93 
Calculated from analysis, corrected for CaCO 
Corrected for CaCO 
Kilocalorie per mole of Ca(OH 
2.41 percent of MgO assumed equivalent to 3.63 percent of CaO 
. : ‘ . . ‘ _ . 
samples can be considered to be formed by the | this investigation. These values are given on the 


equation 


(C/S)\CaO-SiOz- 
H/S)H,.O+AH. (1 


CS)Ca(OH),+SiO,, 0.828H.O0 
H S)H.O+ (C/S+ 0.828 
The sum of the heats of solution of the reactants 
minus the heats of solution of the products is equal 
tothe heat of the reaction. The partial molal heat 
content of the water in HNQOs, 26.3H,O is approxi- 
mately 0.002 keal/mole; in HF, 100H,O it is about 
1.0004 keal/mole [13]. These quantities, which 
represent approximately the heat of solution of the 
Water produced in the reaction, are negligible 
In table 3 are given the calculated heats of reaction 


i Ca(OH), and silica gel to form the samples used in 
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basis of 1 mole of silica. The hydrated calcium 
silicates as occurring in nature or as prepared in the 
laboratory have variable composition. On the 
assumption that the heats of reaction per mole of 
silica would be the same for the substances with 
formulas generally accepted [11], the heats of forma- 
tion of 5CaO-5Si0,-H.O, 2CaO-SiO,-H.O, and 5CaQO.- 
3510,-3H,O from lime, quartz, and liquid water were 
calculated. To perform these calculations, the heat 
of formation [3] of Ca(OH),, SiO,-aq, and liquid water 
were substituted in eq (1) together with the heats of 
reaction as determined from the heat-of-solution 
measurements. The heats of formation of the hy- 
drated calcium silicates were then obtained as the 
sum of the heats of formation of the reactants plus 
the heat of the reaction minus the heat of formation 








of the water appearing on the right side of eq (1). 
Using the heats of formation thus obtained for the 
silicates, thermochemical equations were written for 
the reactions forming 5CaQ-5Si0,-H,O, 2CaQ-siQ.- 
H.O, and 5CaO-3Si0,.-3H.O from lime, quartz, and 
liquid water, and the heats of the reactions cal- 
culated in the usual manner [3]. The values obtained 
are given in table 4. 


TABLE 4 Heats of formation of hudrated calcium silicate 
Bogue Heat of ormat on, A//, 
N € designa Formul Cad 0 
tion |11 mat 
H,OG 
XK onotlite C.SH 5CaO0-5810+-HLO a2 ¢ 
Hillebrandite C,SH(I 2CaO-SiOe-Hy,O 4.7+0.9 
Foshagite Cs°,;H ICaO-3S10,-3H,O 4. ¢ 
The heats of formation from the oxides are calculated he t he se 
ormula The heats of reaction of the reactants Ca(OH ind SiQ»-aq to forn 
the individual compounds was taken from table 3 and the heat rmation of 
the reactants from the oxides calculated from NBS Circular 50 
4. Summary 
The heats of solution of synthetic xonotlite, 


hillebrandite, and natural foshagite have been deter- 
mined in a mixture of nitric and hydrofluoric acids, 
and their heats of formation from the oxides have 
been calculated to 92.6, 34.7, and 94.6 
keal mole, respectively. 


be 
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Effect of Camera Tipping on the Location of the 
Principal Point” 


Francis E. Washer 


Asvmmetric distortion is introduced in a lens-camera combination when the camera is 
incorrectly alined for calibration. A method is described whereby the magnitude and 
direction of the angle of camera tipping can be determined from analysis of the asymmetric 
values of distortion for the case of the object lving at equal but opposite angles from the 


central line of sight. The displacement of the central image from the principal point of 
autocollimation is then readily determined. Theoretical analysis of the problem is given 
and the results confirmed by experiment. The point of symmetry and the principal point of 


autocollimation coincide for the case of no prism effect present on the lens, 


1. Introduction 


One of the problems that arises during the cali- 
bration of precision aerial-mapping cameras is the 
accurate location of the proper principal point [1] * 
with respect to the center of collimation (the point 
determined by the intersection of lines joining 
opposite pairs of collimation index markers). The 
principal point in photogrammetry is defined as the 
point at the foot of the perpendicular drawn from 
the interior perspective center to the plane of the 
photograph. The proper manner of locating this 
point has long been a subject of discussion. There 
are at present three different procedures that are 
employed in locating a point which is believed to 
be a sufficiently close approximation to the true 
principal point. . 

The first method, which has enjoved considerable 
popularity, locates a point called the principal point 
of autocollimation [2, 3], which is the center of the 
image formed in the emulsion plane by the camera 
lens from an incident beam of parallel light, which, 
in the object space, is perpendicular to the emulsion 


plane. This point has usually been designated 


“center cross’”’ [4] by the writer for convenience. 
This point is located by the simplest of the three 
methods, and its popularity undoubtedly derives 
from this simplicity. To locate it, one needs only 
to alm an autocollimating telescope at a distant 
object. The camera ts then interposed between the 
object and the telescope; a plane parallel plate with 
a mirror surface is pressed against the focal plane. 
With the aid of the autocollimating telescope, one 
fan easily adjust the camera until its focal plane is 
normal to the line of sight of the telescope. The 
reflecting plate is then replaced by a photographic 
plate and the distant object is photographed. If the 
collimation index markers are simultaneously photo- 
graphed, one can then easily locate the center cross 
with respect to the center of collimation. 

The second method, which was developed at the 
Bureau, is based on the knowledge that most lenses 


rr 


This is the second paper of a series dealing with problems that relate to the 
bration of precision airplane mapping cameras (see reference [13] at the end of 
8 paper 
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suffer from small errors in centering that cause the 
lens to behave (in a first approximation) as though 
it were composed of a well-centered lens plus a thin 
prism {4, 5]. This is further accentuated in the 
calibration of precision cameras by the presence 
in front of the lens of a filter that is seldom truly 
plane parallel, which accordingly introduces a true 
prism effect. This prism effect causes a displace- 
ment of the center cross from the position it would 
have occupied in the absence of prism effect. This 
immediately poses the question as to which of these 
locations should properly be regarded as the principal 
point. It has been the custom at the Bureau to 
regard the point that would normally have been 
occupied by the center cross in the absence of prism 
effect as the proper principal point of the photograph. 
It is an invariant point when located under standard 
conditions of locating the center cross and is not 
affected by subsequent changes in filters having dif- 
ferent degrees of surface parallelism. The actual 
shift of the principal point from the center cross is 
easily determined from analysis of the asymmetric 
distortion pattern produced by such prism effect. 

The third method came out of attempts to cali- 
brate cameras by a field method [6]. It was not 
customary nor readily practicable to aim the camera 
at the central target of a given target-array with the 
aid of an autocollimating telescope and, conse- 
quently, small errors of alinement sometimes were 
present. Such errors of alinement, called camera 
tipping, can produce an asymmetric distortion pat- 
tern similar but not identical to that produced by 
prism effect. In the analysis of the measurements 
made under these conditions of test, it was found 
possible to make the distortion pattern symmetrical 
by appropriate adjustment of the angles and image 
separations from the central image by computable 
increments. The point about which the distortion 
is tolerably svmmetrical is called the point of sym- 
metry. 

It is evident that for an ideal lens-camera combina- 
tion, Where no tipping of the camera occurs and where 
there is no prism effect in the lens, all three methods 
will lead to the same result. In other words, for 
this case, the principal point of autocollimation (or 
center cross), the principal point, and the point of 








symmetry will be identical. In the present paper, 
the effect of camera tipping on the displacement of 
the central image from the original position of the 
center cross is considered. The magnitude of the 
asymmetric distortion introduced by such tipping is 
evaluated theoretically and confirmed experiment- 
ally. A method of determining the magnitude of 
the angle of camera tipping is shown for the case of 
equal opposite angles. The identity of the principal 
point of autocollimation and the point of symmetry 
is shown for the case of a lens having no prism effect. 


2. Effect of Improper Alinement of Camera 
for Test 


When the camera is being calibrated under condi- 
tions such that the focal plane of the camera is not 
normal to the line drawn to it from the central 
target, the calibration data are likely to show the 
presence of asymmetric distortion [7, 8, 9]. This 
asvmmetric distortion is sometimes confused with 
the asymmetric distortion produced by prism effect 
in the lens [4, 10]. It is, however, much simpler to 
eliminate the asvmmetry of the distortion values 
arising from camera tipping than it is to eliminate 
or reduce the asymmetrical distortion arising from 
prism effect. In this section, the effect of camera 
tipping upon the values of focal leagth and distor- 
tion is considered. 


2.1. Determination of the Equivalent Focal Length 


In figure 1, N is the rear nodal point of an ideal 
lens, L; O is the point where a perpendicular dropped 
from .V intersects the focal plane, consequently VO 





is equal to f, the equivalent focal length. The 
points of intersection with the focal plane of rays 
— 
- o 
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As L is an ideal lens having no 


ON OY f tan B. (] 
Let the camera be rotated by amount ¢ about P 
line through \ normal to the plane of the paper 
The foot of the perpendicular from N moves to 
the image of the central ray is at O’’, and the image. 
is) 


formerly found at X and ¥ are at X’ and Y’, Thy 
relations connecting the various quantities and 
distances are 

OO" =f tan e (9 

O’’X’=f tan(B—e)+f tan e (3 
and 

O”’ Y’=f tan(B+e)—/f tan e. (4 


The distances O’’.X’ and O”’ )” 
ties directly measurable. 
that 


are the only quanti- 
It is clear from figure | 


O”X'+-0"' Y’=f tan (B—e)+f tan(B+e) G 


20 ere meee 2 ftans 
() \ 7 i) } . . . " ’ h 
cos~ e] —tan- B tan?’ e) 
whence 
oO” X’+ i s” : 
= — cos* e¢l—tan-Btan’e). (7 


2 tan B 


For small ¢«, the term (1—&) can be substituted for 
and ¢& for tan’ «. By neglecting the term in 
e‘, eq (7) becomes 


cos” €, 


| oO” xX’ 4 O"Y" 5 3) 
e(i+ tan?’ p)}. Ss 
2 tans 
For very small values of ¢, the term in é& can be 
dropped, and the expression becomes simply 
~ O'X’'+O0"Y'’ 
/ 9 
2 tan Bp 
Thus, the equivalent focal length, f, can be de- 


termined to a good degree of approximation by using 
the measured values of O’”.Y’, O’’Y”’, and the known 
value of 8. The value of f, so determined, will be 
correct to the nearest thousandth of a millimeter fo 
€<10 min for values of the focal length of the orde 
of 150 mm. If eq (9) is used when ¢>10 min, the 
values of f will be higher than the true value, the 
error running as high as 0.012 mm for e=30 min 
In such case it is necessary to determine € by the 
method shown in a later portion of this paper and to 
use this value of €« in eq (8) 
termined with eq (8) agrees with the true value for 
| values of € as high as 30 min 
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2.2. Asymmetric Distortion Produced by Improper 
Alinement 


When a camera is tipped in the manner described 
in section 2.1, it is found that an asymmetric distor- 
tion is introduced by such tipping. This distortion 
affects the relative location in the -image plane even 
gs does the distortion inherent in the lens arising 
from lens aberrations. — There Is, however, one pro- 
nounced difference. Whereas distortion arising from 
lens aberrations is symmetric about the optical axis, 
distortion arising from camera tipping is decidedly 


asymmetric. — : oe . 
In this section, the magnitude of this asymmetric 


distortion will be derived in terms of 6 and e for the 
case of a lens that is distortion-free as far as lens 
aberrations are concerned. By referring to figure 1, 
when distortion is computed from the measurable 
quantities in the usual manner, the distortion, 2);, on 
one side is 


D,=0O" X'—OX (10) 


or 


D), =f|tan(B—e)+ tan e—tan pg}, (11) 


which may be written 
tan e—tan p) 


D ftan B tan e , (12) 
ie l-+tan® tan e 
which for small values of € is equivalent to 
D,=fe tan Ble(1 + tan’ B)—tan 8(1+&)], (13) 
or on neglecting the term in é&, 
fe tan p 
dD fe tan- B+ . (14) 
l ‘ € COs’ 3 
The distortion, //,, for the other side is 
D.=0"' Y’—OY 15) 
or 
D,=f|tan(B+ «)—tan e—tan §}, (16) 
which for small values of « can be show n to be 
; fe“ tan B - 
D.= fe tan’ B+ ——— © (17) 
Cos” 8 


On inspection of the equations for J), and D),, it 
is obvious that for small « the second term containing 
€ will be very small compared to the first term in 
both equations. It is, therefore, clear that D, is 
approximately equal to J), in magnitude but is oppo- 
site in sign. The effect of camera tipping is, accord- 
ingly, an introduction of a symmetric distortion 
whose magnitude is a function of 8 and «. It is 

| Worthy of mention that, when one is analyzing the 
} results of measurement made on a camera whose 


direction of tipping is unknown, the sign of the 
distortion values can serve as a guide in determining 
the angle of tipping and the resultant direction of 
displacement of the central image from the true 
position of the center cross or point of symmetry. 
The point of symmetry or true position of the center 
cross lies on the same side of the central image as 
those points showing maximum negative distortion. 
| To express it another way, the central image is 
| displaced from the center cross toward those points 
showing maximum positive distortion. To illustrate 
the effect of camera tipping, computed values of 
D, and PD), are listed in table 1 for the case of a 
| camera tipped through an angle e=20 min. Even 
| for this relatively large amount of tip, it can be seen 
| that J), is nearly equal to J), but is opposite in sign 
| 





TABLE 1. Asymmetric distortion, D, and Dy, induced in the 
image plane by tipping a camera, equipped with a lens having 
an equivalent focal length of 150 mm, through an angle 


e= 20 min 
Distortion 
Average asymmetric Resultant average distortion 
8 distortion in image 
, 
D D: D.—D Dot Dy 
2 2 
deg mm mm mm mm 
7.5 0.014 0.016 0.015 0. 001 
15 O61 4 063 O02 
22. 5 147 152 . 150 002 
30 287 205 291 004 
37.5 HOT 520 514 006 
5 S63 SAS S73 OLo 


a. Evaluation of the Distortion 


The values of the distortion assigned to a camera- 
lens combination are usually obtained by averaging 
the values measured at the two corresponding 
opposite angular separations from the axis along a 
single diameter. This corresponds to the evalua- 
tion of 


DD, 


we (IS 
9 1S) 


dD 


In the present case, the contribution of lens aberra- 
tions to the distortion has been assumed zero, so if 
D is not zero, the departure therefrom is a conse- 
quence of the camera tipping alone. It is evident 
that for small values of « 


| D,+D, fétang 9 

2 cos* B (19) 

Values of )) have been computed for selected values 

of « and are shown in table 2. It is readily apparent 

that camera tipping does induce a resultant distor- 

| tion that is not eliminated by the averaging process. 

This error in distortion produced by camera tipping 

can be safely neglected for values of ¢ less than 10 

min but may assume serious proportions for values 
of e greater than 20 min. 
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TaAPLE 2. Distortion induced ina normally distortion-free lens 
by tipping the camera, equipped with a lens having an equi- 
valent focal length of 150 mm, through a small angle ¢ 


Distortion for « 
| 
1) min 20 mir jy mu 
deq min mm nim nin 
7 oo oO oot Oo. OO imi 
l ooo mr O08 O13 
22 ml Mme ce 22 
$ ml m4 wu Os 
$4 on (nn ol4 ‘ 
ry O08 on (23 Ol 
For the case of a real lens having measurable 


inherent distortion, it may be inferred that the 
distortion at a given angle 6 is given reasonably 
closely by eq (18) provided that ¢ is kept less than 
15 min. If one is dealing with a lens which is 
required to show less than 0+ .020 mm of distortion 
referred to the equivalent focal length, it is possible 
to induce an error of +0.003 mm in the distortion 
referred to the calibrated focal length if ¢€ is as great 
30 min, It is also likely that small additional 
errors may arise in averaging for those values of 6 
where the distortion characteristics of the lens are 
changing rapidly with 8. Even in such instances, 
the errors in distortion will usually be small com- 
pared with the total distortion 


as 


b. Evaluation of the Angle of Camera Tipping 


When a camera is tested under conditions such that 
the test targets are symmetrically located about the 
central target, it is possible to evaluate e, the angle 


of camera tipping, in a very simple manner. It can 
be shown for small values of ¢ that 
D,—D ' 
“==—Jetan’p 20) 


When )),, 1)),, f, and 8 are known, € can be computed 
from the relation 


(D, dD, (2) 

2f tan* 8 a 
This expression is, of course, computed for the case 
of a distortion-free lens but can be used quite satis- 
factorily for a lens having distortion arising fron lens 
aberrations. This apparent anomaly results from the 
fact that the lens aberration is symmetrical about 
the axis and has the same value at +8 and at —8. 
Consequently, when the difference, )),— 1), is formed, 
the contribution to this difference from lens aberra- 
tion is negligible, the only real contribution being 
the difference in lens distortion between the angles 
8+-e and B—e, which is usually very small for small 
values of «. This difference can become appreciable 
in regions where the distortion is changing rapidly 
with angle. In such cases, a correction can be made, 
if one knows the usual values of the distortion for 
the given type of lens, which usually will be satis- 
factory. Table 1 lists the value of (2),—D,)/2 for 
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one case of camera tipping. If € is evaluated from thp 
values listed, it will be found that the above approxi. 
mation (eq (21)) vields the correct value of «. 

The value of the equivalent focal length, f, useq 
in eq (21) is the approximate value determined With 
the aid of eq (9). The departure of this value of ; 
from its corrected value will usually be too small { 
affect significantly the accuracy of the determination 
of e«. However, the value of ¢, derived from eq (2] 
is sufficiently accurate for use in eq (8) to determine 
the correct value of the equivalent focal length f. 


2.3. Point of Symmetry 


It is clear from the foregoing discussion that if g 
camera is not properly alined for test and if com. 
putations of equivalent focal length and distortion 
are nonetheless performed in the usual manner 
marked asymmetries in the distortion pattern wil 
appear. The magnitude of the asymmetry will, of 
course, depend upon the amount the camera is tipped 
from the condition of proper alinement. In the field 
calibration method when the camera is not initially 
so alined that the focal plane is normal to the line 
connecting the front nodal point of the lens and the 


center target, an appreciable amount of asymmetric 


distortion is introduced. Methods have been de- 
veloped that tend to reduce the magnitude of th 
asvmmetric distortion. This usually involves locat- 
ing a new point in the image plane with respect to 
which the distortion pattern becomes more near 
symmetrical. This point is called the “point of 
symmetry” [6, 11]. 

In the absence of prism effect, point O in figure | 
is the principal point of the camera, and the distance 
0’O"’ is the shift of the central target image from the 
principal point. As O’O’’ =f tan e, and as the valu 
of f tan e can be deduced from the measured value of 
the asvmmetric distortion for a given value of 8 
the location of the principal point can be determined 
For this condition the principal point and point of 
symmetry are identical 


When both prism effect and camera tipping are | 


present, the value of f tan ¢ gives the location of the 
point of symmetry for a selected value of 8. For 
these conditions, the principal point and the point 
of symmetry do not coincide. 


3. Experimental Verification for Tipped 
Camera 


In the calibration of precision cameras by the field 
method, it is usual to locate a point about which the 
distortion is symmetrical. This point is called the 
point of symmetry [6] and this point is sometimes 
set into coincidence with the center of collimation 
instead of setting the principal point into coincidence 
with the center of collimation. The prevailing prac- 
tice at most calibrating laboratories is to set either 
the principal point or the principal point of auto 
collimation [2] into coincidence with the center 
collimation (or fiducial center). The principal pout 
is located at the foot of the perpendicular dropped 
from the center of the exit pupil for the paraxial rays 
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to the image plane and is the point customarily | placed in position on the calibrator. An optical 
brought into coincidence with the center of collima- | plane parallel was set on the focal plane and the 
tion by the Bureau. The principal point of auto- | camera was tipped until its focal plane was normal 
collimation is the image of an infinitely distant point to the optical axis of the autocollimating telescope. 
located in the object space on a line perpendicular | The measured deviation of the line of sight produced 
to the focal plane of the camera. This point is | by the prism was 0.2583°. The tip of approximately 
identical to the point referred to as the center cross | 15 min accordingly may be expected to produce its 
in calibration made by the Bureau. In the absence maximum effect along the diagonal that is the inter- 
of prism effect, the point of symmetry, principal | section of the focal plane and the plane defined by 
point, and the center cross should coincide. The | the optical axes of the collimators in banks I and II. 


truth of the foregoing assertion is readily proved In making the first negative (hereafter referred to 
experimentally and the balance of this section is | as plate 1A), the camera was so oriented in azimuth, 
devoted to a simple experimental proof. that marker A (which designates the line of flight) 
lay between the radial rows of images formed by the 

3.1. Experimental Arrangement camera lens of the targets contained in the collima- 





| tors of banks II and IV. Because of the ti ping, all 
of the collimation markers are moved with respect 
to the target images toward quadrant I. Accord- 
ingly when the location of the central image is deter- 
mined with respect to the center of collimation, O, 
it will be displaced from the center of collimation 
along the diagonal toward the corner of quadrant II. 
The amount of this tip was measured directly and 
found to be approximately 0.2583°. In making the 
second, or check, negative (plate 2B), the camera 
is rotated 180° about the optic axis of the lens with 
all other conditions remaining unchanged. Accord- 
ingly for plate 1A, the central image should appear 
in the quadrant AOD, and displaced from 0 in direc- 
tion II, while for plate 2B, the central image should 
appear in the quadrant BOC and displaced from O 
in the direction II. In the 180° rotation the direc- 
tion IL remains unchanged and AOD is replaced 
by BOC. 


The camera selected for test was one which on 
calibration has been found to have negligible prism 
effect. The separation of the principal point and 
the center of collimation was less than 0.010 mm; 
and the separation of the center cross and the prin- 
cipal point was too small to be measured. This 
camera was placed on the camera calibrator and 
tipped about a horizontal axis lying in a plane defined 
by the optical axes of the collimator banks III and 
IV (see figures 2 and 3 for diagrammatic sketches of 
the trace of this plane in the focal plane of the 
eamera). To achieve a measured amount of camera 
tipping, a small prism was mounted in front of the 
autocollimating telescope, and the telescope alined 
to point at the center of the reticle in the central 
collimator of the camera calibrator. The prism was 
then removed. The camera under study was then 
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FicurE 2. Schematic drawing of test negative obtained on Ficure 3. Displacement of the center image from the center of 
camera calibrator showing the position of the collimation collimation produced by a camera tipping of approximately 
markers A, B, C, and D and the images of the various colli- 15 min. 


mator targets. on : ; , 
_ lhe upper sketches show the displacements computed from the distortion 
The four radial banks are indicated by the numerals I, II, III, and IV. This measurements. The lower sketches show the displacement resolved along the 
figure is to be considered in connection with plate 2B, figure 4. coordinate axes defined by the collimation index markers. 
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3.2. Results of Measurement 


The distances separating all images on the negative 
from the image of the central collimator target were 
measured and are listed in tables 3 and 4 under the 
heading r. (The image of the central collimator 
target corresponds to 0’ in figure 1, and the values 
of r are the measured values of 0’’X’ and 0” Y’ for 
a series of values of 8.) The symbols at the top 
and bottom of the columns headed / tan 8 indicate 
the orientation of the row of images with respect to 
the collimator banks and the fiducial markers in the 
focal plane (see figs. 2 and 3). The values of equiva- 
lent focal length were determined with the aid of 
eq (9) using the measured values of 8, and the value 
of ft for each diameter appears al the top of the 


TARLE 3. 
angle 8 from the center image for plane 1A 


Measured values of distances se parating tmages at 


(deviation of normal by 0.5 prisn 


Camera tipped by approximately 15 min 


diopter D is the distortion for each angle 
f=153.368 mm f=] iu 

( B B D 

8 I II! 
r ftan g D tan 8 D 

lee mm mm mim min nmin en 
4 152. 368 1453. 197 0. 829 153. 249 l 50S 0. 249 
7 117. 234 117. 2 2s 117. 784 117. 78 OO 
0 BAL S1S SA. 464 151 8S, O55 SS. 612 41 
22 t th) (3. 467 O78 63. 58 ( 572 Ol 
15 41. 023 41. O62 asy 41. 12¢ $1. 12¢ TLD 
7 0). 168 20. 184 O16 0, 208 20). 202 OOl 
fh) 0] i ou o { i 
ri 20. 184 20. 168 0. O16 2. 204 1). 20. 0. 001 
15 41.121 $1. O61 (it 41.127 11. 108 oly 
22 63. (33 3. 474 1 63. 609 3, i O46 
0 SS. 775 SS. 467 08 BS. O05 SS. O16 079 
7 118. 029 117. ew 117. 848 17 s O70 
5 153. 690 153. 20 s5 15 7s LA 12 134 

\ > ( 4 

II [\ 
T APLE 4, Measured values of distances se parating images at 


angle 8B from the center image for plate 2B 


Conditions identical with those for 1A except camera is rotated 180° about its 


optical axis 


f=153.341 mm f=153.356 mn 
D \ A ( 
8 I Ii! 
r ftan g D r ftan gs D 
deg mm mm mim mm mm mim 
45 152. 346 153. 170 0. 824 153. 220 153. 54 0. 224 
7. 5 117. 222 117. 542 320 117. 793 117. 782 O11 
10 SS. 200 SS. 449 150 SR. 655 RR. B10 045 
22.5 63. 386 63. 456 070 63. 586 63. 570 O16 
15 41. 020 41. 055 085 41.127 41.125 002 
7. é 20. 164 20. 181 O17 20. 204 “0. 202 002 
0 0 0 0 1] 0 it 
7.5 20. 181 20. 164 0. O17 1). 202 20. 204 0. 002 
15 41.121 41. 054 067 41. 120 41.107 013 
22.5 63. 680 63. 468 167 63. 601 63. 561 040 
10 kA. 700 &S. 452 317 8A. 52 88.614 O68 
37.5 118. 026 117. 543 483 117. 820 117. 775 O54 
iv 153. 701 158. 178 523 153. 356 153. 508 —, 152 
B Cc D B 
Il I\ 








| proper column. By using the proper f for each soy 
| of measurements, the values of f tan 6 are computed 
and the value of the distortion, D, obtained from 
the relation 


D=r—f tan 8. (29) 
These values are shown under the headings f tan g 
and D. It is clear from the measurements that the 
principal effect has been produced along the diree. 
tions | and II as planned. The asymmetric distor. 
tion appearing along the diameter I and II is shown 
graphically as curve 1 in figure 5. It may be noted 
that D, (hereafter referred to as 1) and Dj, (hereafter 
designated II), although opposite in sign, are of 
unequal magnitudes. This difference in magnitude 
is, in large part, a consequence of the large value of 
distortion inherent in the lens itself. 


a. Evaluation of the Angle of Camera Tipping 


Tables 5 and 6 show the evaluation of the angle of 
camera tipping ¢ from the distortion arising from the 
camera tipping. Values of f tan ¢€ are computed 
from the relation 


ttan e221 , 
an € (23 
2tan-B ’) 


derived in section 2.2 (b). For convenience in com- 
putation, the average values of tan 8 and tan? 8 for 
the actual measured angles are listed in table 7. 
The values of f tan ¢ are quite constant for the larger 
angles of 8 for any one set of computations. For 
plate 1A the values range from 0.658 mm at 45° to 
0.691 mm at 30°, with the values of 22.5° and 37.5 
falling within this range. The simple average for 
these four values is 0.678 mm, with the maximum 
departure being —0.020 at 45° for a deviation of 
approximately 3 percent. This spread of values 
can be appreciably reduced if one corrects for the 
presence of distortion, which in this case would 
raise the value of f tan ¢ at 45° by 0.018 to 0.676 
while producing changes not exceeding +0.002 mm 
for any of the other three values. Although this 
would substantially reduce the spread, this correction 
changes the average by less than 1 percent; it does 
not therefore appear worthwhile to include it in the 
present discussion. The values of f tan e at the 
smaller values of 8 are less reliable, as a small error 
ia ¢ produces a large error in f tan e; for example, an 
error of +0.001 mm in r at 8=7.5° produces an 
error in f tan e of +0.029 mm. 
tan ¢€ at 7.5° and 15° accordingly are not used in 
evaluating the average value of f tan «. This type 
of error decreases rapidly with increasing 6 and 
becomes negligible for angles greater than 30°. 
Some of the small discrepancies at the larger angles 
doubtless arise as mentioned from the small differ- 
ences in lens distortion at +e and B—e which are 
small but noticeable for large «. These errors are, 
however, small in comparison to the larger values of 


if tan «. 
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Computation of f tan « or displacement of the center 
induced by camera tipping from distortion values in 


TABLE 5. 
image 
Table 3 for plate 1A 


B I Il Ii—I ftane 
— 
| 
deg mm mm mm mm 
45 0. 829 0. 485 0. 657 0. 658 
37. 5 328 466 397 676 
30 151 308 230 691 
22. 5 —. 078 . 159 118 688 
15 039 060 050 69S 
7.§ 016 016 O16 924 
Average (22.5° to 45°) / tan e=0.678 mm, tan « 
0.004421, and «=0.2533°. | 
8 Ill L\ IV—III ftane 
” 
| 
deg mm mm mm mm | 
45 0). 249 0. 134 0. O58 0. 058 
37 ool O70 036 O61 
+) 041 079 O19 O57 
22.5 o13 046 O16 093 
15 O00 oly OOS 111 
7. 5 ool Ool OOl OSS 
Average (22.5° to 45°) / tan «=0.067 mm, tan e 
0.000437, and «=0.0250 
Resultar tan e=0.681 mm, resultant tan e« 
0.004440, and e=0.2544°. 
- ' 
TaBLE 6. Computation of J tan € or displace ment of the center 


mage induced by camera tipping from the distortion values in 


table 4 for plate IB 


f tan «, 


8 I II li—I tan e 
deg mim mm mm 
{ 0. 824 0. 52 0. 674 0. 67 
37. 5 320 183 . 402 684 
$0) 150 317 234 703 
22 5 O70 167 118 . O8Y 
O35 067 . 051 711 
7 O17 O17 O17 . 982 





} 

} 

| 

Average (22.5° to 4 j/ tan e=0.088 mm, tan e | 

0.004486, and e€=0.2570°. | 

| 

8 Ii! I\ IV—III / tan e 

| 

} 

deg mm mm mm mm | 
45 0. 224 0. 152 0. 036 0. 036 
37.5 Ol 054 . 022 037 
30 045 068 O12 036 
22.5 O16 040 ~O12 070 
15 Oo2 O13 006 O83 
7.5 002 O02 002 115 

Average (22.5° to 45°) / tan e=0.045 mm, tan « 
0.000293, and «=0.0168°. | 
Resultant f tan «=0.689, resultant tan e=0.004492, | 
and e=(.2574°. 


Although care was taken to ensure that the 
maximum effect of camera tipping would take place 
ilong diagonals I and II, the setting in azimuth was 
not quite perfect so a small effect of the tip mani- 
fested itself in directions III and IV. The compu- 
lation of this component of f tan ¢ is shown in the 
wer portions of tables 5 and 6. The full magnitude 
offtan ¢ is then obtained by finding the square root 
ithe sum of the squares of the two orthogonal 
wmponents. The resultant f tan e for plate 1A is 


/'.681, and the corresponding value for plate 2B is 


37 


TABLE 7. Values of average tangent and tangent squared for 
use in interpreting data from camera calibrator 


The measured values of 8 are used in computing tan @. 
— ————4 
l | 
8 tan 8 tan 98 
' 


COLLIMATOR BANKS TAND II 


deg 

7.5 0. 1315530 0. 0173062 
15 . 2677322 . 0716805 
22. 5 . 4138446 . 1712674 
30 . 5768210 . 3327225 
37.5 . 7665412 . 5875854 
45 . 9989116 . 9978244 





COLLIMATOR BANKS III AND IV | 


7.5 0. 1317395 0. 0173553 
15 . 2681092 . 0718825 
22. 5 . 4144998 . 1718101 
30 5778216 . 3338778 
37.5 7680005 . SSGRBR6 
45 1. OOO9S800 1, 0019610 





0.689. It is evident that excellent agreement exists 
between the values, even though the camera was 
rotated 180° and the direction of the effect thereby 
reversed. 

The values of f tan e computed from the measure- 
ments of plates 1A and 2B are the amounts by which 
the center image has been displaced from the position 
it occupies when the camera under test has been 
properly oriented by autocollimation. As the center 
image was initially in near coincidence with the 
center of collimation, the computed displacement, 
can be regarded as a displacement with 
respect to the center of collimation. The magnitudes 
and directions of f tan e¢ derived from tables 5 
and 6 are shown in the upper part of figure 4 for 
plates 1A and 2B. In the lower part of the figure, 
the displacements, originally determined along the 
diagonals, have been resolved into their components 
along AB and CD, which form the reference system 
afforded by lines connecting opposite pairs of col- 
limation index markers. The resultant, 7?, for each 
plate is shown at the bottom. It is clear that there 
is remarkably good agreement between the values 


of R for plates 1A and 2B. 





1A 28 
A 8 
rr mz I | mz 
rl 
0.530 mm 
| 
jones mm 
8) + = GE ——(C o— “T > c 
0.456mm 
ae 
mr I mw I 
a 8 
R=0.699 mm R=0.697mmM 


ACTUAL DEVIATION OF LINE OF SIGHT 0.2563° 
COMPUTED R=0.69I1mm 


FicguRE 4. Measured displacement of the center image from the 
center of collimation on plates 1A and 2B. 


The computed FP is the resultant displacement expected for the measured 


amount of camera tipping introduced. 








Finally, the positions of the displaced center image 
were measured on the plates themselves with respect 
to the center of collimation. The results of these 
measurements are shown in figure 4. Also shown in 
the lower part of figure 4 is the actual measured 
value of the induced angle of camera tipping, which 
is 0.2583°, and the computed value of f tan ¢ derived 
from this value of «. Comparison of these results 
with those shown in figure 3 demonstrates that the 
displacements of the center image computed from 
the asymmetrical values of the distortion agree in 
direction and magnitude with the observed shift 
that actually took place. The small discrepancies 
that remain may be attributed to several factors 
that include the presence of a small residual prism 
in the lens, small additional distortions induced by 
tipping not completely compensated, small amounts 
of plate curvature, and small errors in r and 8. 

Accordingly the resultant f tan « has been de- 
termined in three ways with the following results: 
(1) f tan e¢ has been computed from the measured 
values of the asymmetric distortion, and its average 
value for plates 1A and 2B is 0.685 mm; (2) / tan 
e« has been measured directly on the negative with 
respect to the center of collimation, and its average 
value for plates 1A and 2B is 0.698; (3) finally, f tan 
e has been computed on the basis of the known value 
of the induced camera tipping, and its value is 
0.691. In view of the relatively large value of / tan 
e, it is clear that excellent agreement exists among 
the three methods. This accordingly justifies the 
use of eq (23) in evaluating the angle of camera 
tipping from the asymmetric values of distortion 
and provides a simple method of locating the point 
of symmetry on the basis of distortion measurements. 


b. Evaluation of the Distortion 


In point of symmetry operations, it is customary 
to equalize the distortions. This can be done in the 
present case by subtracting the distortion produced 
by tipping from the measured values. ‘This is done 
in table 8 for one diagonal of plate 1A. In table 8, 
column 1 lists the values of the measured asymmetric 
distortion for one diagonal, and column 2 lists the 
average asymmetric distortion computed for f tan 
«=0.678 mm. By adding the values in 2 to those in 
1, the resultant lens distortion is obtained, which is 


shown in column 3, The distortion values in 3 are 
nearly equal for corresponding values of 8. This 
process is shown graphically in figure 5. In the 


figure, curve 1 shows the values of the distortion 
computed with the displaced central image as the 
center. Curve 2 shows the amount of this distortion 
attributable to camera tipping. Curve 3 shows the 
amount of distortion still remaining after the cor- 
rection for camera tipping is applied. It is evident 
that curve 3 is already quite symmetrical and closely 
approximates a typical distortion curve for this 
variety of wide-angle lens. In the table, column 4 
shows the complete equalization obtained by aver- 
aging. Column 5 shows the difference between 3 
and 4 and leads to the belief that complete equal- 
ization was not achieved because of the introduction 
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TABLE 8. Adjustment of distortion values from tippe 


data for plate 1A, collimator banks I and J] 


EFL=153.368 mm, / tan «=0.678 mm, CFL 


d camer, 


153.310 mm 


8 l 2 3 ] t) ° 
deg mm mm mm mm mm mm mm 
45 —{(), 829 0. 678 —(. 153 —(). 172 0.019 0.058 |~9 14 
37.5 -. 328 3YS . 070 069 ool 044 13 
30 —. 151 226 . 075 078 —. 003 033 lll 
22. 5 —. 078 116 . 038 040 O02 024 064 
15 O39 . 049 O10 O10 . 000 016 . 028 
7.5 016 O12 —. 004 000 004 OO8 008 
0 0 0 0 0 0 0 0 
7.5 O16 —(). O12 004 000 O04 
15 O60 —. (49 oll O10 ool 
22.5 159 116 043 . 040 003 
30 308 226 OS2 O78 004 
37 i} SUS . 068 069 - OO1 
45 185 676 191 172 O19 


* Column 

Initial asymmetric values of distortion. 

f tan e tan ? 8 (contribution from camera tipping 

Value of distortion compensated for camera tipping 

Averaged values of distortion referred to equivalent focal length 
5. Departure of compensated distortion from average 

6. Aftan 8, where Af=0.058 mm. 

7. Values of distortion referred to calibrated focal length 


of a small amount of residual distortion induced by 
plate tipping. Column 6 shows the correction jy 
distortion introduced by changing from the EFL to 
the CFL by amount Af=0.058 mm. Column 7 
shows the values of the distortion referred to the 
calibrated focal length. These final values of the 
distortion referred to the CFL are nearly the same gs 
those obtained from the properly oriented camers: 
the small differences arising from the uncompensated 
distortion introduced by camera tipping, which in 
this instance would not exceed +0.003 mm. 
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Ficure 5. Evaluation of the true values of distortion from the 


asymmetric measured values. 


The open circles (curve 1) show the actual values of distortion as initially 
determined from the measurements. The triangles (curve 2) show the valueé 
asymmetric distortion produced by a tipping of amount «=0.2533° (see table § 
or 15.2min. The closed circles (curve 3) show the symmetrical distortion patt® 
remaining when curve 2 is subtracted from curve 1. 
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4. Discussion 


When a camera under calibration is initially so 
aiined that the focal plane of the camera is not 
normal to the line joining the front nodal point of 
the lens and the distant target, the values of the 
distortion based on the assumption that the central 
image is indeed in its proper position will be asym- 
metrical about this point. The values of the dis- 
tortion at points located at equal but opposite 
angles from the center can be used to evaluate the 
magnitude and direction of the displacement of the 
central image from its position for the condition of 
true normality. On making this correction, which 
can be done in a simple manner, the point of sym- 
metry can be readily located. For the case of no 
prism effect on the lens, the point so located is the 
point of autocollimation or center cross. Accordingly 
it is clear that for these conditions, the point of 
autocollimation and point of symmetry are identical 
and can be regarded as the true principal point of 
the camera. 


The author expresses his appreciation to other 
members of the staff for assistance rendered during 
the course of this work. The negatives used in the 
confirmatory experiments were made by W. P. 
Tayman. Measurements on the negatives were 
made by W. R. Darling. The drawings were made 
by E. C. Watts. 
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Mass Spectra of Some Lead Alkyls 
Edith I. Quinn, Vernon H. Dibeler, and Fred L. Mohler 


Monoisotopic mass spectra are 
methylethyllead, dimethyldiethyllead, 
trimethyl-n-butyllead, 
butyllead. 


given 
methyltriethyllead, 
trimethyl-sec-butyllead, 
Relations between relative probabilities of ion formation and molecular structure 


for 


the tri- 
tetrapropyllead, 


and __ tetra-iso- 


compounds: tetramethyllead, 
tetraethyllead, 


trimethyl-tert-butyllead, 


are observed for several ion configurations and are discussed briefly. 


1. Introduction 


The systematic comparison of the mass spectra . 
a group of related compounds not only affords : 
basis for chemical analysis but also provides infor- 
mation on the molecular structure and _ relative 
probabilities of bond dissociation processes. Organ- 
ometallic compounds are particularly interesting for 
this purpose, and studies of mass spectra have been 


| 
| 


reported for a number of these, including the tetra- | 


methyl compounds of the group IV elements (1), 
lead and mercury alkyls [2], and the eyclopentadieny] 
compounds of V, Cr, Fe, Co, Ni, Re, Ru, Mn, and 
Mg [3]. 

The concomitant data of isotope abundance ratios 
obtained for the metallic elements in organometallic 
compounds is also of interest. Although possible 
interference from ions formed by rearrangement 
during the ionization-dissociation process (e. g., 
hydride ions) preclude their use for discovering rare 
isotopes, satisfactory agreement has been obtained 
in the comparison of routine analyses of lead, using 
the halides and the tetramethyl compound [4], and 
of mercury, using the element and the dimethyl 
compound [5]. 

The present study, made possible through a gift of 
the lead alkyls from the Ethyl Corp., Detroit, 
extends the published mass spectral data of lead 
alkyls to include trimethylethyllead, dimethyldi- 
ethyllead, woes beat lar ad, tetra-n-propyllead, 
trimethyl-n-butyllead, trimethyl-sec-butyllead, — tri- 
methyl-fert-butyllead, and tetra-iso-butyllead. 
ative dissociation probabilities are compared for the 
several molecules and conclusions noted with regard 
to molecular structure and mass spectra. 


2. Experimental Details 


Mass spectra were obtained by using a Consoli- 
dated 21-103 mass spectrometer and conventional 
procedures. In preliminary experiments, liquid 
samples were introduced through a gallium-covered 
porous glass frit [6]. However, considerable holdup 
was Observed for the heavier liquids, and as a result 
the sample pressure did not reach a steady value and 
pump-out time was excessive. The use of a mercury- 
sealed orifice [7] avoids these difficulties. There was 
nd indication of reaction between the mercury and 
the lead alkyls. 


Figures in brackets indicate*the"liter ature"references at the’end“of his paper 
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A resolving power of about 1 in 500 was obtained 
by reducing the width of the collector slit from 32 to 
6 mils and by the use of a high potential on the 
“metastable suppressor.”’ This electrode is located 
between the collector slit and the collector, and in 
addition to filtering out ions that have lost energy 
through metastable transitions, it also acts as a 
virtual limiting slit. 

Table 1 lists the lead alkyls and the purities as 
stated by the supplier. The tetramethyllead and 
tetraethyllead were part of the same samples used in 

previous study [2]. Recognized impurities in the 
samples were principally alkyl halides, except in the 
case of the trimethyl-tert-butyllead, which contained 
tert-butanol in excess of the total reported impurity 
and may have resulted from some decomposition of 
the sample. Mass spectral contributions of these 
impurities were subtracted from the lead alkyl 
spectra. Impurities consisting of other organo-lead 
compounds are not easily identified. One check on 
the purity is to compare the spectrum of the vapor 
in contact with an excess of liquid with the spectrum 
obtained when a small sample of liquid is completely 
vaporized. In the first spectrum, lower vapor- 
pressure impurities are suppressed whereas those 
with higher vapor pressures are enhanced. This 
technique disclosed approximately 0.1 percent of 
trimethylethyllead and dimethyldiethyllead in the 
tetramethyllead sample. 








TaBLE 1. Purity of lead alkyls as reported by the supplier 
— — — 
Compound Formula Purity 

Trimethylethyllead | (CH )3PbC.Hs 99* 
Dimethyldiethyllead (CHs3)2Pb(C2Hs)2-. Wy 
Methyltriethyllead | CHsPb(C2Hs5)s- 99r 
Tetra-n-propyliead (n-C3H7)4Pb_~ | 95 
Trimethyl-n-butlyllead (CHs3)sPb(n-CyHp) - -- 95 
Trimethyl-sec-butyllead (C H)3Pb(sec-C4He) -- 95 
Trimethy]-ert-butyllead (CH )3Pb(tert-C Hy) 99 
Tetra-iso-butyllead (i-C4Ho)«Pb 95 

The relative abundance of lead isotopes in 

natural ores is subject to large variations. ‘“Com- 


mon lead,’’ which is frequently a mixture from sev- 
eral sources, may also exhibit variations. In this 
research, the isotope ratios are computed inde- 
pendently for each compound. However, as the 
preparation of the Prone dates to 1934 and the 








origin of the leads is unknown, small variations in 
the isotope ratio observed in these samples, although 
probably real, are not a matter of scientific interest. 


3. Results 


The polyisotopic spectra of the lead alkyls have 
been submitted for publication elsewhere [8]. As 
dissociation processes are difficult to interpret from 
mass spectra of molecules containing elements with 
a number of isotopes of comparable abundance, 
monoisotopic spectra are computed. First, contri- 
butions of ions containing C’ and D are subtracted 
in the conventional manner by using a C'*/C" ratio of 
91.0 and an H/D ratio of 6700. Common lead is com- 
posed of four isotopes with the ate percent- 
age abundances: 204:206:207:208=1.3:26:21:52. 
Thus each ion of formula Pb(R), gives rise to four 
mass peaks in addition to the peaks resulting from 
ions containing C" and D atoms. <A method for ob- 
taining the isotope abundance ratios and applying 
them to the computation of the monoisotopic spectra 
has been described in detail [1, 

Table 2 summarizes the prince ipal ions in the 70-v 
monoisotopic spectra of 10 lead alkyls, including 
tetramethyl and tetraethyllead. The abundances 
are relative to the most abundant ion taken as 100. 
Observed molecule ions are in boldfaced type. The 
sensitivities (ion current per unit of pressure in the 








sample reservoir) of the most abundant ions are ex- | 
pressed as relative to the sensitivity of the C,H? | 
ion of n-butane. These appear in the last line of 
the table. 


4. Discussion 


The hydrocarbon ions in general, and partic ularly 
the me thy] ions in the lead alkyl mass spectra, may 
exhibit variations in relative abundance for instru- 
mental reasons, and conclusions based thereon are 
subject to some uncertainty. Nevertheless, some in- 
teresting results are observed. The relative abun- 
dance of the CH} ions in tetramethyllead, the meth. 
ylethylleads, and tetraethyllead consistently de- 
creases with the number of Pb-CH, linkages. [yn 
the methylbutylleads the CH} abundance is af.- 
fected by the configuration of the butyl radical as a 
low CH? abundance occurs for the n-butyl, whereas 
a high abundance is observed for the tertiary butyl 
compound. 

The abundance of C,H;* ions also follows a logical 
relation to molecular structure, except in the case of 
tetraethyllead. Formation of the C,H;* ion by 
dissociation of the Pb—C,H, linkage is relatively 
less probable than formation of CH,* from Pb(C H;),. 
~The relation of abundance of propyl and butyl ions 
to molecular structure is evidently complicated. 
The formation of C,H;* ions is relatively improbable 
in the tetrapropy] compound, and analogous to the 
tetraethyllead. ‘The increased abundance in the 
trimethyl-tert-butyllead is unexpected because of 
the complicated dissociation process involved and 
the low abundance in the tetra-iso-butyllead. The 
abundance of the C,H,* ion in the sec-butyl com- 
pound is not consistent with the order of bond 


energies: primary >secondary > tertiary. 


TABLE 2. Principal ions in the monoisotopic spectra of some lead alkyls * 
Ion Pb(R!)4 Pb(R R?2) | Pb(R')2(R2)»| Pb(R!)(R? Pb(R?), Pb(R), Pb(R R*) Pb(R R Pb(R'),(Ré Pb(R’), 
CH 5. 35 3. 44 2. 06 1.42 0. 21 1. 45 1.02 2. 85 12.5 3. 58 
CoHs - 0. 96 1. 50 1. OS 1. 57 1. 82 4.78 5. 57 0.31 2. 21 
C 3H: 1. 16 0. 47 1.41 3.19 0. 85 
C.He * 97 10. 2 2.M4 5.17 
Pb 69. 7 76. 1 78. 1 72.6 46.2 85.3 76.4 67.4 67.1 86. 8 
PbH 3. 24 16. 6 0.1 13.5 39. 2 32.8 4.3 21.2 74 $3. 7 
PbCH; 6. 35 5. 93 4.44 1. 3€ 0. 63 6.42 6.12 4.68 ). 38 6.00 
PbCH 87.8 100.0 95. 1 78. 6 00 100.0 9.1 87.1 87. 5 100. 0 
PbC 2H, 0.72 1.11 27.5 0.23 0.16 1.1 
PbC,H 0.53 11.8 31.3 60. 100.0 14.7 61 79 0.91 41.0 
PoC. « 25.0 15.0 3. 58 0.14 2.79 13.8 10.4 4.0 11.9 
PbC;H 3.12 4. 4 
PbC3H 0. 27 0. 38 1.43 1. 63 0. 26 14.0 0.3 3. OF 
PbC\H SI 10.0 If 17.9 9, 69 2. 22 13 
PbC gH, 100. 0 81.3 4.12 1.42 62.0 100.0 100.0 100.0 49.3 
PbC«Hs 0.19 0.19 1.04 36. 5 9. 24 0. 28 1. 61 
PbC,H, 1. 67 0 14 3. 4 24 
PbC.Hy, 44.4 100.0 11.9 3. 24 34.7 0. 32 0. 1¢ 1 
PbCsHy 1s 100.0 0.43 23.3 67 15 39 
PbCsHy, 13 0. 66 , 
PoC Hy 0.09 76.7 10.0 68.0 7.0 0. 36 
POHCcH ys 1, 06 10.4 0.18 0.27 0. 30 
PoC: Hy 3. 02 1. 21 
PoC:Hiw 0. 96 0. 46 0.17 1.79 0, 10 0. 38 
PbCysHy 7.31 45 
PbCyHa 24 
Sensitivity relative 
to n-butane m/e= 
43 0. 69 0.60 0.89 0. 86 1.2 0. 72 0. 76 0. 94 0.9 0. 89 
* R'=CHs, R*’e@C,;Hs, R'=C3H:, R'=n-CyHy, Ri=s-C,H», R®=t-CyH», R?=iso-CyHy 


42 





art 
ha 


res 
ful 


of 
alt 
rel 
me 


of 

the 
mo 
ra- 
the 
spe 


col 
the 
rad 
low 
pot 
tet 
are 
ti01 
fig 


are 
Ine 
rad 
abu 
and 
Pbt 
leac 


abu 
The 


in t 





arly 
nay 
tru- 
are 
 in- 
un- 
*th- 
de- 

In 
af- 
iS a 
reas 
ityl 


ical 
p of 
by 
ely 
[;),. 
ons 
ted, 
ible 
the 
the 

of 
and 
The 
ym- 
ond 





Hydrocarbon ions apparently resulting from re- 
arrangement during ionization-dissociation processes 
have been noted previously for tetramethyllead [1]. 
Traces of impurities in the compounds used in this 
research effectively masked suc h ions and prevent 
further discussion concerning their origin. 

Lead atom ions are of comparable abundance in all 
of the compounds. The abundance of PbH* ions, 
although comparable for most of the compounds, is 
relatively low in the tetramethyllead and the tri- 
methyl-tert-butyllead. 

The PbCH,* ion is one of the principal ions in all 
of the lead alkyl spectra except tetraethyllead. Al- 
though the latter requires rearrangement to form the 
monomethyllead ion, the tetrapropyllead and tet- 
ra-iso-butylleads also require rearrangement, and 
the PbCH,* ion is the most abundant ion in the 
spectra of the ‘se compounds. 

The abundances of the PbC,H,* ion are generally 
consistent with the structure of the lead alkyls, e. g., 
the abundance increases with substitution of ethyl 
radicals for methyl groups of tetramethyllead. The 
low abundance of this ion in the methylbutyl com- 
pounds is consistent with structure. However, the 
tetrapropyllead and _ tetra-iso-butyllead compounds 
are exceptions. Formation of this ion by dissocia- 
tion of a hydrogen atom from a dimethyllead con- 
figuration evidently occurs with low probability. 

The abundances of the dimethyllead ion, PbC,H,*, 
are also generally consistent with molecular structure. 
Increased substitution of ethyl radicals for methyl 
radicals in the methyl-ethyl leads results in decreased 
abundance. Formation of this ion in tetrapropyllead 
and the butvl leads is by rearrangement. The 
PbC,H,* ions are evidently hydrides of the dimethyl- 
lead ion formed by rearrangement. 

As expected, the trimethyllead ion, PbC,H,*, 
abundant for molecules having 3 or 4 methyl groups. 
The relatively large abundance of this ion observed 





in the tetrapropvllead and tetra-iso-butyllead spectra | 
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dihydride of the 
This is suggested by the apparent 
high stability of the ion configuration, PbR,*, in 
which three bonding orbitals of the central lead 
atom are occupied by hydrogen atoms or alkyl 


suggests a monopropyllead 


H,PbC,H,* 


groups. Such stability is demonstrated by the large 
abundances of the Pb(CH;),*+ ions of tetramethyl- 
lead, trimethylethyllead, and the trimethylbutyl- 
leads, the Pb(CH;).C.H;+ ion of dimethyldiethyl- 
lead, ‘the Pb(CH;,) (CH; )o + ion of methyltriethyllead 
and the Pb(C,H;),;* ion of tetraethyllead. Tetra- 
propyllead and _ tetra-iso-butyllead are evidently 
exceptions to this, probably for steric reasons. How- 
ever, the trimethyl-n-butyllead spectrum includes 
an abundant ion, PbC,H,,;*, which very probably has 
the trialkyl configuration, Pb(CH;).C,Hyg. 

lons heavier than PbC,H,,;* are relatively rare. 
Molecule ions of only five of the lead alkyls were 
observed with abundances as large as 1 percent of 
the most abundant ion. This fact complicates the 
mass spectrometer analysis of mixtures of these 
compounds, particularly for those containing small 
amounts of lower molecular weight compounds. 
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Frequency Response of Second-Order Systems 
With Combined Coulomb and Viscous Damping’ 


Thomas A. Perls? and Emile S. Sherrard 


1. Introduction 


The increasing importance of vibration measure- 
ments for both military and nonmilitary applica- 
tions has stimulated the production of a large variety 
of vibration-measuring instruments in recent years. 
This activity in design and manufacture has not 
been matched by progress in the analysis of the 
response Of vibration instruments, or in the exten- 
sion or modernization of previous analyses of their 
response. It is still frequently assumed, in the 
application of these instruments, that the response 
of an actual vibration instrument is identical with 
the response of an ideal instrument in which the 
damping is entirely ‘‘viscous,”’ i. e., proportional to 
the velocity of the mass or “‘seismic element.”? The 
actual response may differ significantly from the 
ideal response if the damping varies with other than 
the first power of the velocity, or is dependent on 
displacement, or if any appreciable amount of 
coulomb damping is present. This last type of 
damping, sometimes referred to as dry or sliding 
friction, exerts on the moving element a force that 
is constant in magnitude and always acts in the 
direction opposite to the velocity of the moving 
element. 

Coulomb damping is usually present in a mechani- 
eal second-order system such as is illustrated in 
figure 1. Here, in addition to the coulomb damping, 
a spring force and a viscous damping force are shown 
as acting upon the seismic element. The motion of 
such a system has been analyzed by Den Hartog * 
in & 1931 paper in which results are presented 
graphically in a series of figures. For application to 
vibration instruments in present use, Den Hartog’s 

This work was conducted under a program of basic instrumentation research 
ud development sponsored by the National Bureau of Standards, the Office of 


Naval Research, the Air Research and Development Command and the Atomic 
Energy Commission. 


om” at the Missile Systems Division ,Lockheed Aircraft Corp., Van Nuys, 
alif. 

: J, P. Den Hartog, Forced vibrations with combined coulomb and viscou 
metion, Trans. Am. Soe. Mech. Eng. 53, APM 53 9, 107-115 (1931). 














Curves obtained with an analog computer are presented for the magnification factor 
versus frequency ratio of second-order systems with combined coulomb and viscous damping. 
The ranges of the parameters are as follows: Viscous damping ratio from 0.05 to 5.0, in 15 
steps; coulomb damping ratio from 0 to 0.9, in 11 steps; and frequency ratio from 0 to 2.0. 
Boundaries between regions with 0, 1, and 2 stops per half-cycle are also shown. 


results require extension to larger values of viscous 
friction and a larger range of frequencies. 

The present paper extends these results through 
the ranges applicable to such inertial instruments as 
accelerometers and jerkmeters. These inertial in- 
struments are usually designed for viscous damping 
between 0.6 and 0.7 of critical damping, and are 
usually employed between zero frequency and a 
maximum frequency less than the resonance fre- 
quency (of the mass and spring). Particular empha- 
sis is therefore placed in this paper upon results for 
viscous damping between 0.6 and 0.7 of critical 
damping, and for frequencies between zero and the 
resonance frequency. 

Results are presented graphically in figures 2 to 
17. These curves extend from zero to twice resonance 
frequency, and are presented so as to be readily 
applicable to vibration instruments used below their 
resonance frequency. They were determined without 
particular difficulty through the use of an analog 
computer. The original program called for the 
computation of an additional set of curves for 
frequencies between twice the resonance frequency 
and 50 times the resonance frequency. These 
curves would have been applicable to devices such 
as seismometers and velocity meters, which are 
ordinarily employed at frequencies above their 
resonance frequency. Unfortunately, the analog 
computer accuracy proved poor at these high fre- 
quencies. The computation of displacement at high 
trequencies consisted essentially of the double inte- 
gration of the difference of a sinusoidal forcing func- 
tion and the coulomb friction. When the analog 
computer attempted this double integration, its results 
contained appreciable errors caused by drifting and 
hunting. The size of these errors caused a termina- 
tion of the high-frequency computations. Hence, 
only curves applicable to vibration instruments 
usually employed below their resonance frequency 
are presented in this paper. 


45 








2. Differential Equation for Vibration 
Instruments 


figure 1 shows a typical vibration instrument 
with its frame rigidly attached to a sinusoidally 
vibrating structure having a motion X cos wf. The 
displacement, x, of the seismic element relative to the 
frame of the instrument satisfies the differential 
equation: 
(1) 


Miz+Cir+ Keri F=MXe’ cos oat, 


where (see fig. 1) 


M= mass of seismic element. 

(’=viscous damping force per unit velocity. 

K=restoring force exerted by the spring per 
unit displacement. 

F=coulomb damping force, assumed constant 

in magnitude and changing sign so as 

always to oppose the motion. 

X=amplitude of forcing motion. 

w= 22f=angular frequency of forcing motion. 


A number of equations in Den Hartog’s paper give 
a complete literal solution of this equation with the 
right-hand side of the equation replaced by a term 
P cos (wt+¢), where P is the amplitude of a ‘‘periodic 
disturbing force.” For practical application Den 
Hartog expressed his results by a set of computed 
curves of magnification factor versus frequency ratio 
for seven values of viscous damping ratio between 0 
and 0.5 of critical damping, and for a number of 
values between 0 and 0.9 of the ratio of coulomb 
damping force to peak disturbing force. Most of 
Den Hartog’s curves extend over a frequency range 
of 0.5 to 2.0; the one for zero viscous damping ex- 
tends over a frequency range of 0 to 2.0. Curves for 
zero coulomb damping (i. e., viscous damping only) 
are also available elsewhere.‘ In the present paper, 
all the curves extend over the frequency ratio range 
between 0 and 2.0. They are given for 15 values of 
viscous damping ratio between 0.05 and 5.0, and for 
11 values of coulomb damping ratio between 0 and 
0.9. 


3. Summary of Results 


The curves shown in figures 2 to 17 are drawn 
through points plotted by an analog computer. Each 
figure presents a family of curves for a constant value 
of C/C., the ratio of viscous damping (, to critical 
viscous damping, C,=2VKM. Each 
family is a plot of magnification factor or acceleration 


curve of the 
. . 2; rarele , " / 

response Zow;/X versus the frequency ratio, w/w,, 

where 2 is the maximum value of the displacement z, 


w,=yK/M is the undamped resonance frequency of 


‘See, for example, C. 8. Draper, W. McKay, and S. Lees, Instrument engineer- 
ing. Il, Mathematics (McGraw-Hill Book Co., Inc., New York, N. Y., 1953 
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the mass-spring system, and X=w*r is the peak 
acceleration of the forcing motion. For inertia) 
instruments, the term “acceleration response” seems 
appropriate for row,?/X because this quantity may be 
interpreted as the product of two important instry. 
ment characteristics: 29/X, 
ment) of the seismic element per unit acceleration 
applied to the instrument, and w;, the resonance 
frequency of the instrument, to which instrument 
sensitivity is inversely proportional. 

Each curve of the family is drawn for a particular 


the response (displace. 


value of the parameter ay/X, where ay is the minij- 
mum acceleration that must be applied to the instrp- 
ment to overcome coulomb friction and _ produce 


motion. This parameter ap/X is equivalent to the 
parameter F/P used by Den Hartog. Either may be 
obtained from the other by employing the relations 


F= Ma, and P=MX. 

The dashed line plotted in each figure is a boundary 
line between steady-state motions with no stop and 
one stop per half-cycle. In the region above this 
boundary, the steady-state motion is continuous 
without stops. In the region immediately below this 
boundary, the motion stops once each half-cycle. 

Only these two regions, nonstop and one-stop 
motion, were computed for C/C,>0.2 and w/w, >0.1, 
For C/C,<0.1, four regions characterized by no stop, 
one stop, two stops, or more than two stops each 
half-cycle were computed. When more than two 
regions occur, the magnification factor versus fre- 
quency ratio curves become very irregular for /a, 
less than about 0.4. For w/w, <0.1, larger and larger 
amounts of computing time would be required for 
accurately determining the magnification factor as 


C/C. becomes smaller and smaller. Hence for 
C'C'.=0.2, 0.1, and 0.05, the curves are terminated 
at w/w,—0.1. 


4. Accuracy of Results 


The curves of figures 2 to 17 are plotted to permit 
reading the magnification factor to within 1 percent 
of the full scale of each figure. These curves were 
drawn through points plotted by an analog computer. 
Over-all accuracy of the points plotted by the analog 
computer, which consisted of elements with an ac 
curacy of +0.1 percent of rated output, is estimated 
at +0.3 percent of the full-scale magnification factor. 
A number of check points computed manually and 
by the computer agreed with each other to within 
+1.0 percent. Total computer errors and drafting 
errors are believed to introduce, in any figure, imac 
curacies of no more than 1 percent of the full-scale 
magnification factor of each figure. 

It is believed that greater accuracy in comput 
tion would have resulted from the use of a digital 
computer that employed Den Hartog’s literal equa- 
tions. However, an estimate of total costs showed 
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analog computation to be more economical. One 
reason was the availability in the analog computer’s 
auxiliary equipment of a plotting board, which could 
partially process the results of the computer study by 
jotting points of the magnification-factor curves. 
The other important important advantage of the 
analog computer was its use of eq (1) rather than the 
involved literal equations of Den Hartog’s paper, 
which would have been employed by the digital com- 
puter. Use of the analog computer avoided the ex- 
penditures necessary for a careful scrutiny of these 
equations for possible misprints. 


5. General Effects of Coulomb Friction Upon 
Instrument Response 


The presence of coulomb friction makes instru- 
ment response dependent on the amplitude of the 
sinusoidal excitation. This amplitude dependence 
isshown graphically in figure 18, which is a cross-plot 


of the data for C/C,.=0.65, as given in figure 10. In 
figure 18, magnification factor is plotted versus 


X/ay for various constant values of w/w,. 

The asymptotes shown as horizontal lines in the 
right-hand portion of figure 18 represent the magni- 
fication factors at vibration amplitudes sufficiently 
large to make negligible the effect of coulomb friction. 
An inspection of the figure shows that, for C/C,=0.65, 
X/ay must be >25 if the actual response is to be 
within 5 percent of the asymptote for all values of 
WS Wp. 

Rules similar to this can be stated for other values 
of C/C.. For C/C.>0.3 and for 0<w/#,<1.0, a 
general rule of thumb for estimating the error intro- 
duced by coulomb friction is to assume that this 
error is equal to the error at w/w,=0. According to 
this general rule, a given percentage value of a,/X 
results in that same percentage reduction of magnifi- 
eation factor from its ideal (a¢/X=0) value. This 
rule is a good approximation for the reduction in re- 
sponse for C/C, in the range 0.6 to 0.7. As the value 
of C/C. moves further and further outside this range, 
the approximation becomes poorer and_ poorer. 
However, the damping in most inertial instruments 
is between 0.6 and 0.7 of critical damping; hence, 
this rule is useful for practical inertial instruments. 

A physical example of the effect of coulomb frac- 
tion upon instrument response is shown in figure 19. 
This figure reproduces an oscillogram from two 
different velocity meters mounted back to back and 
subjected to the same sinusoidal forcing function. 
The upper trace shows that one instrument is stop- 
ping once each half-cycle as the result of coulomb 
friction. The lower trace shows the other to be 
exhibiting essentially ideal behavior. 

Figure 20 shows the time variation of displacement 
and velocity for various values of ap/X, for a low 
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value of viscous friction (C/C.=0.1), and a low value 
of frequency ratio (w/w,=0.1). These records were 
made during the computer study. Each set of two 


traces is for a different value of ap/X. At ay/X=0.3, 
one-stop motion occurs, as evidenced by the long dwell 
of the velocity trace on the (central) zero line. For 


ay/X equal to 0.5, the second trough of the damped 
sinusoidal oscillation in the velocity trace approaches 
the zero line. At ap/X=0.525, this trough reaches 
the zero, and the motion stops twice during each 
half-cycle. Further increase of ap/X to 0.55 fur- 
nishes an additional restraint, which prevents 
another start after the motion stops at nearly the 
same point of the velocity trough. Thus we have 
again one stop per half-cycle. The motion remains 
one-stop until displacement and velocity decrease 
sufficiently for a stop to occur at the first trough of 
the velocity, giving two-stop motion once more. 
This occurs at approximately ay/X=0.8. However, 
an increase of ap/X to 0.9 furnishes the additional 
restraint necessary to prevent another start. Thus 
motion again becomes one-stop. 

As C/C, decreases from 0.1, the damped sinusoid of 
velocity becomes less damped and it is conceivable 
that for w/w, <0.1, motion may stop not only at the 
third, second, or first trough of the velocity, but also 
at a fourth, fifth, or higher-order trough of the 
velocity. In such cases, the response of the instru- 
ment becomes very irregular. 


a 


The analog computer study was performed by W. 
McCool and B. Zimmerman of the Naval Research 
Laboratory. The authors thank Messrs. McCool 
and Zimmerman for their suggestions and excellent 
cooperation in this study. M. Abramowitz and R. 
Dressler supplied information on the use of the digital 
computer. W. A. Wildhack continually furnished 
stimulating ideas during the investigation. A. Bucek 
prepared the curves and illustrations shown in this 
report. 
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Figure 7. Instrument response for C/C, 
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Fiaure 8. Instrument response for C/C, 
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Figure 9. Instrument response for C/C,.=0.60. 
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Figure 10. Instrument response for be 
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Ficure 11. Instrument response for C/C,=0.70. 
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FicurE 12. Instrument response for C/C, 
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Figure 13. Instrument response for C/C.=0.80. 
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Ficure 14. Jnstrument response for C/C.=0.90 
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Figure 15. Instrument response for C/C.=1.0. 
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Ficure 16. Instrument response for C C= 8.0. 
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Figure 17. Instrument response for C/C.=5.0. 
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FIGURE 18. Effect of amplitude upon instrument response, 
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FIGURE 19. Ex pe rimental response of two different 
mbration instruments. 
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